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ABSTRACT 

We analyze the Type II Plateau supernovae (SN II-P) 2005cs and 2006bp with the non-LTE model 
atmosphere code CMFGEN. We fit 13 spectra in the first month for SN 2005cs and 18 for SN 2006bp. 
Swift ultraviolet photometry and ground-based optical photometry calibrate each spectrum. Our anal- 
ysis shows both objects were discovered less than 3 days after they exploded, making these the earliest 
SN II-P spectra ever studied. They reveal broad and very weak lines from highly-ionized fast ejecta 
with an extremely steep density profile. We identify Hen4686A emission in the SN 2006bp ejecta. 
Days later, the spectra resemble the prototypical Type II-P SN 1999em, which had a supergiant-like 
photospheric composition. Despite the association of SN 2005cs with possible X-ray emission, the 
emergent UV and optical light comes from the photosphere, not from circumstellar emission. 

We surmise that the very steep density fall-off we infer at early times may be a fossil of the com- 
bined actions of the shock wave passage and radiation driving at shock breakout. Based on tailored 
CMFGEN models, the direct-fitting technique and the Expanding Photosphere Method both yield 
distances and explosion times that agree within a few percent. We derive a distance to NGC 5194, the 
host of SN 2005cs, of 8.9±0.5Mpc and 17.5±0.8Mpc for SN 2006bp in NGC 3953. The luminosity of 
SN 2006bp is 1.5 times that of SN 1999em, and 6 times that of SN 2005cs. Reliable distances to Type 
II-P supernovae that do not depend on a small range in luminosity provide an independent route to 
the Hubble Constant and improved constraints on other cosmological parameters. 

Subject headings: radiative transfer - stars: atmospheres - stars: supernovae: individual: 2005cs, 
2006bp, 1999gi - stars: distances - stars: evolution 



1. INTRODUCTION 

Although radiation represents a negligible fraction of 
the energy in core-collapse supernova (SN) explosions, 
most of what we know about these events is inferred 
from the analysis of spectra and light curves. The grav- 
itational binding-energy of the newly-formed protoneu- 
tron star is on the order of 10 53 erg. Within millisec- 
onds after core bounce, a few x 10 52 erg of this energy 
is used to photodissociate the infalling outer iron core. 
Starting with an electron-neutrino burst when the core 
reaches nuclear densities, the radiation of neutrinos of 
all flavors operates over a few tens of seconds after the 
bounce as the protoneutron star cools, and carries away 
a few x 10 52 erg. These neutrinos are believed to deposit 
~10 51 erg (1% of the total) into the infalling progenitor 
mantle, reversing the accretion, and to provide the in- 
ternal and kinetic energy for the SN ejecta (Woosley & 
Janka 2005). Only -10 49 erg (0.01% of the total) is even- 
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tually processed into light, that gets radiated with a rate 
equivalent to a few x 10 8 L© sustained for three months. 

Depending on the progenitor radius, the shock emerges 
at the surface between a few hours and a day after core 
bounce, with a radiative precursor that is expected to 
heat and accelerate the surface layers, producing a soft 
X-ray flash (Chevalier 1982; Ensman & Burrows 1992; 
Matzner & McKee 1999; Blinnikov et al. 2000). Cool- 
ing, due to adiabatic expansion and radiative losses, is 
moderated as time progresses by recombination of hydro- 
gen during the photospheric phase of Type II SN and by 
non-thermal excitation by unstable isotopes. Although 
UV emission is significant for the first few days after 
shock breakout, the spectral energy distribution (SED) 
subsequently peaks further and further into the red, first 
in the visual and then in the near-IR, a few months af- 
ter explosion (Kirshner et al. 1973,1975; Mitchell et al. 
2002; Leonard et al. 2002a; Brown et al. 2007). 

Core-collapse SN explosions constitute an excellent 
laboratory for inferring the properties of massive stars at 
the end of their lives and, indirectly, their pre-SN evolu- 
tion. As the material expands, the photosphere recedes 
to deeper and deeper layers of the star's former enve- 
lope, allowing the observer to probe all the way from the 
surface at shock breakout to the innermost layers just 
outside the compact remnant. The nebular phase, the 
epoch when the ejecta become optically thin throughout 
(~3-4 months after explosion for Type II SNe), permits 
the inspection of the regions where the blast originated, 
offering a means to constrain the mechanism and the 
morphology of the explosion as well as the nucleosyn- 
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thetic yields. 

Detailed quantitative analyzes focus on the early, pho- 
tosphcric, evolution (Hoeflich 1988; Eastman & Kirsh- 
ner 1989; Schmutz et al. 1990; Baron et al 1996, 2000, 
2003, 2004, 2007; Mitchell et al. 2002; Dessart & Hillier 
2005a, 2006ab), because the ejecta are best suited for 
"standard" model atmosphere computations: There is 
an optically-thick base where the radiation thermalizes; 
there is no contribution from non-thermal excitation by 
isotopes like 56 Ni and its daughter product 56 Co; electron 
scattering dominates the opacity; and line and contin- 
uum formation is spatially confined. Advances in com- 
puter technology permit better handling of non-Local- 
Thermodynamic-Equilibrium (non-LTE) effects and the 
critical role played by metals is now systematically ac- 
counted for (Baron et al. 2004; Dessart & Hillier 2005a). 
Non-LTE effects are important even at early times, for 
example, to reproduce the Hei lines observed in optical 
spectra. 

Here, we present a quantitative spectroscopic analysis 
of the early photospheric-phase evolution of the Type II 
SNe 2005cs in NGC 5194 and SN 2006bp in NGC 3953, 
both with plateau-type light curves as shown in Fig. [TJ 
We employ the non-LTE model atmosphere code CMF- 
GEN (Hillier & Miller 1998; Dessart & Hillier 2005a), 
in its steady-state configuration, following a similar ap- 
proach to Dessart & Hillier (2006a; hereafter DH06) and 
their analysis of SN 1999em. A novel component of our 
study is the use of Swift (Gehrels et al. 2004) Ultra- 
Violet/Optical Telescope (UVOT; Roming et al. 2005) 
ultraviolet photometry that places important constraints 
on the SED in the UV, the early cooling of the photo- 
sphere, and the reddening. We infer from our models 
that SN 2005cs and SN 2006bp, were discovered only -2 
days after explosion, the earliest detections after explo- 
sion since SN 1987A (Type II peculiar) and SN 1993J 
(Type lib; see, e.g., Matheson et al. 2000). As a result 
of this prompt discovery, we have the opportunity to ob- 
serve a bright UV spectrum and to see high-excitation 
lines that have not previously been observed in SN II-P. 

Finally, the fine time sampling of our spectra during 
the first month, along with dense photometric monitoring 
allows us to make an accurate distance determination. 
This is a good demonstration of the potential of Type 
II-P SNe, when coupled with detailed models, to provide 
accurate extragalactic distances. Unlike the standard ap- 
proach that depend on the Cepheid period-luminosity re- 
lation in the LMC and other nearby galaxies, our method 
is based on physical models fitted directly to observa- 
tions, and the distances it provides are independent of all 
the lower rungs of the distance ladder. We aim to assess 
the accuracy of distance determinations based on Type 
II-P SNe, and to refine the value of the Hubble constant 
Ho, in the spirit of the earlier efforts of Schmidt et al. 
(1994). Measurements of the total matter and baryonic 
content of the Universe via CMB experiments or large- 
scale structure correlations depend on the square of Hq 
(Spergel et al. 2007; Eisenstein et al. 2005). In concor- 
dance solutions, these uncertainties affect measurements 
of the Dark Energy density and its associated equation- 
of-state parameter w (Garnavich et al. 1998; Riess et 
al 2005; see Macri et al. 2006 for an illustration of the 
effect an improved accuracy on Hq has on w). A small 
improvement in our knowledge of Hq will lead to a sig- 



nificant improvement in the understanding the energy 
content of the Universe. 

This paper is structured as follows. In the next section, 
we present the photometric and spectroscopic datasets. 
In [JH we present the model atmosphere code CMFGEN 
and describe our methods. In $4j we present a quan- 
titative spectroscopic analysis at 13 photospheric-phase 
epochs for SN 2005cs, and in Sj5j analyze 18 epochs for 
SN 2006bp. In fj6j we apply the Expanding Photosphere 
Method (EPM) as well as a direct-fitting method (in 
the spirit of the Spectral-fitting Expanding Atmosphere 
Method, SEAM, of Baron et al. 2000) to infer the dis- 
tance and the explosion date for both SN events. In £0 
we discuss the implications of our results, and in fJSl we 
conclude. 

2. SUMMARY OF OBSERVATIONS AND DATA 
REDUCTION 

To sample the early photospheric phase with as many 
optical observations as possible, we use multiple data 
sources. For both objects, Swift UVOT photometry is 
used. Spectra are less frequently observed than the pho- 
tometry, typically once every few days, and this sets the 
limit to our analysis. In the next sections, we review 
the ground-based optical photometry, the optical spec- 
troscopy, the Swift UVOT data. 

2.1. Optical photometry 

UBVr'i' photometry was obtained for SN 2005cs and 
SN 2006bp at the FLWO 1.2m telescope on Mt. Hop- 
kins. The Minicam instrument was used for SN 2005cs 
while the Keplercam was used for SN 2006bp (for details 
on cameras and photometric reduction see Hicken, 2007 
in preparation). All data were bias subtracted and flat 
fielded using standard procedures. PSF-fitting photome- 
try was performed using Dophot (Schechter et al. 1993). 
Host-galaxy subtraction was performed for SN 2006bp. 
The host-subtracted photometry was approximately 0.02 
mag fainter than the photometry from the unsubtracted 
images, across all bands. A host-galaxy subtraction tem- 
plate was not available for SN 2005cs but the supernova 
was on a relatively faint and smooth background and 
observations were taken while SN 2005cs was bright so 
the absence of host subtraction is negligible. Compari- 
son stars for each SN were calibrated on only one night. 
However, our confidence in their accuracy is high as other 
SNe that had been observed on other photometric nights 
were observed on these nights, and the results are con- 
sistent. The comparison stars of these other SNe showed 
good agreement across all nights (the standard deviation 
of a typical comparison star across these nights being 
about 0.015 mag in V). We present the photometry ob- 
tained for SN 2005cs in Tableland in Fig.[T](left panel), 
and for SN 2006bp in Table [2] and in Fig. CD (right panel), 
where, for completeness, we have also included, where 
appropriate, the observations obtained by Pastorello et 
al. (2006; hereafter, P06) and Tsvetkov et al. (2006), 
and the Swift observations for the three bluest filters. 

2.2. Optical spectroscopy 

Some of the spectra for SN 2005cs presented and an- 
alyzed in this paper have been previously published by 
P06, while a subset for SN 2006bp have been presented 
by Quimby et al. (2007). We refer the reader to those 
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Fig. 1. — Observed light curves for SN 2005cs (left; optical data include the observations of T06, P06, and those of the CfA) and for 
SN 2006bp (right). Besides optical magnitudes in UBVRI (i' is shown instead of I for SN 2006bp), we also include the Swift UVOT 
magnitudes in the filters uvw2, uwm2, and uvwl (corresponding values and photometric errors are given in Table [T]& [4] for SN 2005cs and 
Tabl e [2] fe [5] f or SN 2006bp, for each filter and date). [See the electronic edition of the Journal for a color version of this figure, and see 
£|2.1l and i|2.3l for discussion.] 



TABLE 1 

Journal of photometric observations of SN 2005cs 



UT Date 

(- 


JD 

h2,453,000) 




U 


B 






V 




R 






1 






r 






i' 




Source 


2005-06-30 


552.4 


13 


48 


± 0.05 


14.36 


± 





05 


14.48 


± 


0.02 


14.46 ± 





04 


14.44 


± 


0.04 












P06 


2005-07-01 


553.2 








14.42 


± 





02 


14.42 


± 


0.02 














14.68 


± 0.03 


14 


34 


± 0.03 


CfA 


2005-07-01 


553.3 


13 


53 


± 0.06 


14.36 


± 





05 


14.46 


± 


0.05 


14.40 ± 





02 


14.34 


± 


0.06 












P06 


2005-07-02 


554.2 


13 


50 


± 0.02 


14.47 


± 





01 


14.41 


± 


0.01 














14.64 


± 0.02 


14 


32 


± 0.02 


CfA 


2005-07-02 


554.5 


13 


69 


± 0.10 


14.45 


± 





04 


14.51 


± 


0.03 


14.47 ± 





03 


14.29 


± 


0.10 












P06 


2005-07-03 


555.2 








14.52 


± 





02 


14.46 


± 


0.02 


















14 


35 


± 0.02 


CfA 


2005-07-03 


555.3 








14.45 


± 





02 


14.45 


± 


0.02 


14.30 ± 





02 


14.37 


± 


0.05 












T06 


2005-07-04 


556.2 








14.55 


± 





01 


14.52 


± 


0.02 














14.66 


± 0.02 


14 


36 


± 0.02 


CfA 


2005-07-05 


557.4 








14.52 


± 





04 


14.54 


± 


0.04 


14.42 ± 





05 


14.29 


± 


0.08 












P06 


2005-07-06 


557.8 


14 


02 


± 0.06 


14.60 


± 





05 


14.55 


± 


0.03 


14.36 ± 





03 


14.40 


± 


0.03 












P06 


2005-07-06 


558.2 


13 


87 


± 0.02 


14.61 


± 





02 


14.52 


± 


0.02 














14.60 


± 0.03 


14 


33 


± 0.03 


CfA 


2005-07-06 


558.3 








14.57 


± 





02 


14.54 


± 


0.02 


14.34 ± 





02 


14.29 


± 


0.03 












T06 


2005-07-07 


559.2 








14.58 


± 





03 


14.53 


± 


0.01 














14.59 


± 0.02 


14 


33 


± 0.01 


CfA 


2005-07-07 


559.4 








14.65 


± 





04 


14.56 


± 


0.09 


14.41 ± 





09 


14.35 


± 


0.11 












P06 


2005-07-08 


560.2 
















14.54 


± 


0.01 


















14 


23 


± 0.01 


CfA 


2005-07-08 


560.3 








14.72 


± 





01 






























CfA 


2005-07-09 


561.3 








14.74 


± 





02 


14.50 


± 


0.02 


14.26 ± 





02 


14.25 


± 


0.03 












T06 


2005-07-10 


562.2 


11 


64 


± 0.03 


14.84 


± 





02 


14.52 


± 


0.02 
























CfA 


2005-07-11 


563.1 








14.93 


± 





01 


14.50 


± 


0.02 
























CfA 


2005-07-11 


563.3 
















14.50 


± 


0.02 
























T06 


2005-07-11 


563.4 


14 


84 


± 0.05 


14.90 


± 





04 


14.54 


± 


0.04 


14.37 ± 





02 


14.24 


± 


0.05 












P06 


2005-07-12 


564.3 








14.92 


± 





02 


14.50 


± 


0.02 


14.32 ± 





02 


14.20 


± 


0.03 












T06 


2005-07-13 


565.4 








15.08 


± 





05 


14.65 


± 


0.04 


14.42 ± 





04 


14.29 


± 


0.06 












T06 


2005-07-13 


565.4 


15 


27 


± 0.06 


15.04 


± 





05 


14.64 


± 


0.02 


14.37 ± 





03 


14.29 


± 


0.05 












P06 


2005-07-14 


566.3 


15 


36 


± 0.06 


15.19 


± 





02 


14.63 


± 


0.02 


14.39 ± 





02 


14.32 


± 


0.03 












T06 


2005-07-14 


566.4 


15 


29 


± 0.05 


15.10 


± 





05 


14.67 


± 


0.02 


14.42 ± 





04 


14.27 


± 


0.02 












P06 


2005-07-15 


567.3 








15.13 


± 





03 


14.60 


± 


0.02 


14.34 ± 





02 


14.23 


± 


0.03 












T06 


2005-07-17 


569.4 


15 


92 


± 0.07 


15.31 


± 





05 


14.67 


± 


0.03 


14.40 ± 





03 


14.27 


± 


0.04 












P06 


2005-07-19 


571.3 








15.40 


± 





02 


14.72 


± 


0.02 


14.43 ± 





02 


14.27 


± 


0.03 












T06 


2005-07-19 


571.4 








15.34 


± 





02 


14.72 


± 


0.02 


14.38 ± 





04 


14.22 


± 


0.04 












T06 


2005-07-19 


571.4 








15.36 


± 





12 


14.72 


± 


0.05 


14.45 ± 





09 


14.26 


± 


0.07 












P06 


2005-07-20 


572.4 








15.39 


± 





07 


14.71 


± 


0.03 


14.45 ± 





03 


14.27 


± 


0.03 












P06 


2005-07-25 


577.4 








15.46 


± 





09 


14.72 


± 


0.04 


14.36 ± 





06 


14.24 


± 


0.06 












P06 


2005-07-27 


579.4 








15.60 


± 





07 


14.73 


± 


0.04 


14.39 ± 





04 


14.16 


± 


0.04 












P06 


2005-07-28 


580.3 








15.69 


± 





05 


14.73 


± 


0.03 


14.42 ± 





03 


14.20 


± 


0.04 












T06 



NOTE. — Sources for photometric observations arc CfA, P06 and Tsvctkov et al. (2006; T06). For the distance determinations performed in j|6. 11 wc 
use the photometry corresponding to the dates on which spectroscopic observations were obtained. If there is no photometry on the date of a spectrum, 
wc linearly interpolate the slowly varying magnitude obtained at bracketing dates. 
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TABLE 2 

Journal of photometric observations of SN 2006bp 





UT Date 


— 

JD 

( i o a cr q r\c\r\\ 

(+2,453,000) 




- — 

U 






- — 

B 








— 

V 






— 

r 










— 

i 




2006-04-10 


836.27 








15.59 


± 


0.02 


15. 


■14 


± 


0.02 


15.42 


± 





.02 


15. 


.39 


± 


0.02 


2006-04-11 


837.20 


14.97 


± 


0.03 


15.47 


± 


0.02 


15. 


31 


± 


0.02 


15.24 


± 





.02 


15. 


.23 


± 


0.02 


2006-04-13 


839.20 


14.92 


± 


0.03 


15.44 


± 


0.02 


15. 


16 


± 


0.02 


15.03 


± 





.02 


1 1. 


.9 1 


± 


0.02 


2006-04-14 


840.23 








15.43 


± 


0.02 


15. 


12 


± 


0.02 


14.97 


± 





.02 


14. 


.91 


± 


0.01 


2006-04-17 


843.24 


15.04 


± 


0.03 


15.46 


± 


0.02 


15. 


11 


± 


0.02 


14.94 


± 


0. 


.02 


14. 


.87 


± 


0.02 


2006-04-18 


844.16 


15.09 


± 


0.03 


15.47 


± 


0.02 


15. 


12 


± 


0.02 


14.93 


± 


0. 


.02 


14. 


.86 


± 


0.02 


2006-04-20 


846.17 


15.20 


± 


0.02 


15.52 


± 


0.02 


15. 


15 


± 


0.02 


14.93 


± 





.02 


14. 


.87 


± 


0.02 


2006-04-22 


848.22 


15.33 


± 


0.03 


15.59 


± 


0.02 


15. 


18 


± 


0.02 


14.94 


± 





.02 


14. 


.90 


± 


0.02 


2006-04-23 


849.23 


15.42 


± 


0.03 


15.63 


± 


0.02 


15. 


19 


± 


0.02 










14. 


.91 


± 


0.02 


2006-05-02 


858.27 


16.44 


± 


0.03 


16.01 


± 


0.02 


15. 


22 


± 


0.02 


14.93 


± 





.02 










2006-05-03 


859.19 


16.55 


± 


0.04 


16.07 


± 


0.02 


15. 


24 


± 


0.02 


14.93 


± 





.02 


14. 


.89 


± 


0.02 


2006-05-05 


861.17 


16.74 


± 


0.04 


16.16 


± 


0.02 


15. 


27 


± 


0.02 


14.95 


± 





.02 


14. 


.88 


± 


0.02 


2006-05-06 


862.16 


16.79 


± 


0.05 


16.21 


± 


0.02 


15. 


26 


± 


0.02 


14.96 


± 


0. 


.02 


14. 


.88 


± 


0.02 


2006-05-11 


867.19 


17.12 


± 


0.06 


16.39 


± 


0.03 


15. 


32 


± 


0.02 


15.00 


± 





.02 


14. 


.89 


± 


0.02 


2006-05-14 


870.20 








16.48 


± 


0.04 


15. 


37 


± 


0.02 


15.05 


± 


0. 


.02 


14. 


.92 


± 


0.02 


2006-05-30 


886.23 


18.36 


± 


0.07 


16.85 


± 


0.02 


15. 


52 


± 


0.02 


15.11 


± 


0. 


.02 


14. 


.98 


± 


0.02 



Note. — Log of CfA photometric observations. For the distance determination done in S|6.2| we only use the £?, V, and i' magnitudes. We 
use photometric data on the dates of the spectra, or, if there is no photometry on the day of a spectrum, linearly interpolate between magnitudes 
obtained at bracketing dates. 



TABLE 3 

Journal of spectroscopic observations of SNe 2005cs and 2006bp with FLWO 1.5m+FAST 



UT Date a 


JD a 


Epoch b 


APA C 


Airmass 


Flux stds. d 


Seeing 


Exp. time 


Observer^ 




(+2,453,000) 


(d) 


(°) 






(") 


00 




SN 2005cs 


2005-06-30.23 


551.73 


4 


1.1 


1.31 


BD28, BD17 




1200 


RH 


2005-07-01.27 


552.77 


5 


5.0 


1.53 


BD28, BD17 


4 


900 


RH 


2005-07-02.23 


553.73 


6 


3.1 


1.34 


BD28, BD26 


2 


1200 


RH 


2005-07-03.15 


554.65 


7 


32.7 


1.09 


F34, BD33 


2 


720 


RH 


2005-07-04.168 


555.66 


8 


85.4 


1.11 


BD28,BD17 


1 


1200 


JG 


2005-07-04.18 g 


555.68 


8 


78.4 


1.15 


BD28,BD17 


1 


1200 


JG 


2005-07-05.18 


556.68 


9 


8.9 


1.16 


F34, BD17 


1 


1200 


JG 


2005-07-06.17 


557.67 


10 


5.2 


1.15 


BD28, BD17 


1 


1200 


JG 


2005-07-07.15 


558.65 


12 


1.8 


1.11 


BD28, BD17 


1 


1200 


EF 


2005-07-09.19 


560.69 


13 


0.7 


1.24 


BD28, BD17 


1-2 


1200 


PB 


2005-07-10.18 


561.68 


14 


1.1 


1.22 


BD28, BD17 


1-2 


1200 


PB 


2005-07-11.18 


562.68 


15 


2.9 


1.23 


BD28, BD17 


1-2 


1200 


MC 


2005-07-12.16 


563.66 


16 


2.8 


1.17 


BD28, BD17 


1-2 


1200 


MC 


2005-07-28.17 


579.67 


32 


1.0 


1.38 


BD28, BD17 


1-2 


1200 


PB 


2005-07-29.15 


580.65 


33 


2.1 


1.30 


BD28, BD17 


2-3 


960 


MC 


2005-12-25.56 


730.04 


182 


5.0 


1.11 


BD28, BD17 


1 


1800 


MC 


SN 2006bp 


2006-04-20.30 


845.80 


12 


12.3 


1.13 


F34, HD84 


1.5 


1200 


WP 


2006-04-21.21 


846.71 


13 


3.5 


1.07 


F34, HD84 


1.5 


1200 


WP 


2006-04-24.19 


849.69 


16 


67.9 


1.08 


F34, BD26 


1.5 


1500 


JH 


2006-05-01.35 


856.85 


23 


0.8 


1.41 


F34, HD84 


1.5 


780 


WB 


2006-05-02.14 


857.64 


21 


3.8 


1.10 


F34, HD84 


1.5-2 


1200 


WB 


2006-05-03.22 


858.72 


25 


64.4 


1.08 


F34, HD84 


1.5-2 


1200 


TG 


2006-05-04.15 


859.65 


26 


2.8 


1.09 


F34, HD84 


1.5-2 


1200 


TG 


2006-05-07.19 


862.69 


29 


12.2 


1.07 


F34, BD17 


1.5-2 


1200 


WP 



Note. — All the above spectra were obtained with the 1.5 m Tillinghast telescope using the FAST spectrograph (Fabricant ct 
al. 1998). The grating used yields a resolution of ~ 7A. The observed wavelength range roughly spans 3500-7400 A. (See gUJfor 
discussion). 

a UT Date and Julian Date at midpoint of observations.' 3 Epoch of spectrum relative to the estimated explosion datc. c Absolute 
difference between the observed position angle and average parallactic angle over the course of the observation.^ Standard star pairs 
used for flux calibration: BD17 = BD+17°4708; BD26 = BD+26°2606; BD28 = BD+28°4211; BD33 = BD+33°2642; F34 = Feige 
34; HD84 — HD 84937. e Seeing is based upon estimates by the observers/ Observers: EF — E. Falco; JG — J. Gallagher; JH — 
J. Hernandez; MC = M. Calkins; PB = P. Bcrlind; RH = R. Hutchins; TG = T. Groncr; WB = W. Brown; WP = W. Peters. 5 Both 
spectra of SN 2005cs on UT 2005-07-04 were combined into a single spectrum for plotting and modeling purposes. 
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papers for more information on those spectra. In the 
analysis section of SN 2005cs, we indicate this source as 
P06, wherever appropriate, while in the analysis section 
of SN 2006bp, we indicate this source as Quimby. 

All the other spectra were obtained with the 1.5 m Till- 
inghast telescope at the Fred Lawrence Whipple Obser- 
vatory (FLWO) using the low-dispersion FAST spectro- 
graph (Fabricant et al. 1998). Observational details of 
the spectra are listed in Table [3] The FAST data were all 
reduced in the same manner, using standard IRAF 8 rou- 
tines for CCD processing and spectrum extraction. All 
the spectra were extracted using the optimal algorithm 
of Home (1986). The spectra are then wavelength cali- 
brated using calibration-lamp spectra; the calibration is 
later checked against night sky lines and adjusted ac- 
cordingly. We used our own IDL procedures for flux 
calibration. These include the removal of telluric lines 
using the well-exposed continua of the spectrophotomet- 
ric standard stars (Wade & Home 1988; Matheson et al. 
2000). 

The relative spectrophotometry is good to about 5%, 
when done at the parallactic angle (Matheson et al. 2007, 
in prep.), but no attempt was made to determine an 
absolute flux calibration. To avoid second-order light 
contamination, we calibrated each spectrum using two 
standard stars of different colors (see Table [3|) . In gen- 
eral, the spectra were observed at the parallactic angle 
to minimize the effects of atmospheric differential refrac- 
tion (Filippenko 1982). When this is not the case, the 
spectra can suffer from a relative flux depletion at bluer 
wavelengths, especially when observed at high airmass 
(secz > 1.1, where z is the zenith angle). In Table [3] we 
list the angle in degrees, APA, between the observed po- 
sition angle and the parallactic angle for each spectrum, 
along with the airmass. When APA > 10°, the relative 
spectrophotometry is less accurate. 

Around maximum light, both SNe 2005cs and 2006bp 
were bright with respect to their host galaxy. Never- 
theless, the earliest (June 30 th UT) and latest (Dec 25 th 
UT) FAST spectra of SN 2005cs have substantial galaxy 
background, and the relative flux calibration is affected 
by this contamination. 

Finally, all SN 2005cs spectra presented here have 
been de-redshifted, adopting our own measurement of 
466±35km s~ 1 of the recession velocity of NGC5194, 
while for SN 2006bp and NGC 3953, we adopt 
1280 km s _1 (Quimby et al. 2007; Verheijen & Sancisi 
2001; a comparable value of 1053 km s _1 is given by Ho 
et al. 1997). 

2.3. Swift UVOT photometry 

Swift UVOT photometry allows us to track the ultra- 
violet color evolution of SNe from a few days to a few 
weeks after shock breakout. Observations of both SN 
2005cs and SN 2006bp were obtained with Swift UVOT 
(Gehrels et al. 2004; Roming et al. 2005; Poole et al. 
2007) at multiple epochs covering from just hours after 
discovery until the UV light had faded to the background 
emission level of the SN environment. The Swift UVOT 

8 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, operated by the Association of Universities for Research 
in Astronomy, Inc., under contract to the National Science Foun- 
dation of the United States. 



data for SNe 2005cs and 2006bp, originally presented in 
Brown et al. (2007) and Immler et al. (2007), have been 
reanalyzed using the improved calibration of Poole et al. 
(2007). The magnitudes in Table [4] and in Tableware 
in the UVOT system, distinguished by lower case letters. 
The UVOT u, b, and v are close to the Johnson UBV; 
however, the UVOT u band is bluer than the Johnson 
U so there can be significant differences at early times 
when the SNe are blue. The b and v agree well with their 
Johnson counterparts. The magnitudes and flux densi- 
ties were measured using a 5" aperture source region, cor- 
responding to that used to calibrate the zeropoints and 
coincidence loss corrections, using uvotmaghist. Here, 
we subtract the background light by means of late-time 
template images, selecting a background region with the 
same average count rate per pixel as in the 5" source re- 
gion. When the background subtracted count rate drops 
below that of the background (as measured in the tem- 
plate images) we report the magnitude of the background 
as an upper limit. The errors reported include the sys- 
tematic errors in the zeropoints. 

We use the flux equivalent of the Swift UVOT mag- 
nitudes given in Table |4] & [3 accounting for the dif- 
ferent (and varying) energy distribution of the super- 
nova by adding an additional systematic error of 10% 
in quadrature with the statistical error. We plot these 
fluxes at the effective wavelengths of the UVOT filters 
along with appropriately reddened synthetic flux distri- 
butions computed with the non-LTE model atmosphere 
code CMFGEN (see §0 §01 and §©. Swift UVOT fil- 
ters uvw2, uvm2, uvwl, u, b, and v have their effective 
wavelength at 2030A, 223lA, 2634A, 350lA, 4329A, and 
5402 A (Poole et al. 2007). 

3. THEORETICAL METHODOLOGY 

The analysis of SNe 2005cs and SN 2006bp is anal- 
ogous to that presented in DH06, which analyzed SN 
1999em at 8 epochs in the first 45 days after explosion. 
That paper determined its distance by means of various 
EPM-based approaches (Kirshner & Kwan 1974; East- 
man & Kirshner 1989; Schmidt et al. 1994; Eastman 
et al. 1996; Hamuy et al. 2001; Leonard et al. 2002a; 
Elmhamdi et al. 2003; Dessart & Hillier 2005b) and an- 
other method akin to the SEAM (Baron et al. 2000). 
Here, we have two advantages: Swift observations moni- 
tor the UV light and we have more densely sampled data. 
Swift provides additional constraints on the reddening, 
the temperature, and the ionization state of the ejecta 
at early epochs. Our analysis covers 13 epochs in (^30 
days) for SN 2005cs and 18 epochs in (~45 days) for 
SN 2006bp, providing better tracking of the spectral and 
photometric evolution of each. 

Time-dependent terms in the statistical- and radiative- 
equilibrium equations become important a few weeks af- 
ter explosion (Utrobin & Chugai 2005; Dessart & Hillier 
2007a, b) . We avoid that complication by focusing on the 
early-time photospheric phase. We know that assuming 
a steady state in the statistical and radiative equilibrium 
equations produces a prediction of the ejecta ionization 
state that is somewhat below its true level, that Hei 
lines persist over a longer period than we predict, and 
that strong Balmer lines at late times are a clear sign of 
the epoch when this assumption is no longer valid. 

CMFGEN solves the radiative transfer equation in the 
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TABLE 4 

Swift UVOT Photometry of SN 2005cs 



UT Date 


JD 






uvwl 


uvm2 


uvwl 


u 






b 










V 






onnc; c\a on 

2UUo-Uo-oU 


O A CQCCO 
2400002. 


.4 


13 


0/ ± U.03 


i n a A _L n hq 
12.64 ± U.U3 


12. /8 ± O.Oo 














1 A 

14. 


.63 


± 





.03 


onnc; cv? no 
zUUo-U /-Uo 


2453555. 


i 




1 _l_ n HQ 

21 ± U.03 


















1 -1 
14. 


63 


i 

± 


u. 


.Ub 


2005-07-05 


2453557. 


3 


15 


.05 ± 0.04 


14.87 ± 0.04 


14.39 ± 0.04 


13.64 ± 0.03 


14 


.75 


± 





.03 


14. 


68 


± 





03 


2U05-0/-06 


245350 / . 


.9 


15 


26 ± U.04 


lo.lo ± U.U4 


i a cn J— n c\a 
14. bU ± 0.04 


i q 7t; _i_ r\ nQ 
13.75 ± U.U3 


14 


76 


± 


0. 


.03 


14. 


.66 


i 

± 


0. 


03 


2005-07-07 


2453558. 


.8 




>15.7 


>15.7 


14.98 ± 0.04 


13.85 ± 0.03 


Id 


.76 


± 





.02 


14 


,64 


± 


0. 


.03 


2005-07-08 


2453560. 


.0 




>15.7 


>15.7 


15.63 ± 0.05 


14.14 ± 0.03 


14 


.78 


± 





.03 


14 


,59 


± 





.03 


2005-07-09 


2453560. 


6 




>15.7 


>15.7 


>15.7 


14.32 ± 0.04 


Id 


,78 


± 





.04 


14 


,58 


± 





.04 


iuuo-u r-iu 


240o00Z. 







>15.7 


>15.7 


>15.7 


1/1 C.1 _1_ C\ f\A 

14.0/ ± U.U4 


14. 


,87 


± 





,05 


14. 


.63 


± 


0. 


.05 


2005-07-11 


2453562. 


.8 




>15.7 


>15.7 


>15.7 


14.92 ± 0.05 


Id 


,98 


± 





,05 


14. 


.70 


± 





.05 






5 




>15.7 


>15.7 


>15.7 














14. 


.67 




o 


05 


2005-07-14 


2453565. 


.5 




>15.7 


>15.7 


>15.7 


15.43 ± 0.05 


15 


,21 


± 





.04 


14 


.74 


± 





.04 


2005-07-17 


2453568. 


.5 




>15.7 


>15.7 


>15.7 


>15.8 


15 


.48 


± 





,05 


14 


.76 


± 





.05 


2005-07-22 


2453573. 


.5 




>15.7 


>15.7 


>15.7 


>15.8 


15 


,62 


± 





,06 


14. 


90 


± 





.05 


2005-09-14 


2453628. 


5 




>15.7 


>15.7 


>15.7 


>15.8 


16 


.33 


± 


o 


.06 


14. 


95 


± 





.04 












TABLE 5 
































Swift UVOT Photometry of SN 2006bp 






















UT Date 


JD 






uvw2 


uvm2 


uvwl 


u 






b 










V 






2UUO-U4-1U 


z4oooob. 





14 


49 ± 0.03 


14.74 ± 0.03 


14.62 ± 0.03 


14.61 ± 0.02 


15. 


75 


1 

zb 





02 


lb 


63 


1 





03 


zUUo-U4-ll 


2453836. 


.6 


14. 


.57 ± 0.04 


14.83 ± 0.05 


14.66 ± 0.04 


14.60 ± 0.04 


15 


,60 


=t 





.04 


15. 


.52 


± 





05 


2UUO-U4-11 


2453836. 


9 


14. 


.62 ± 0.03 


14.83 ± 0.04 


14.63 ± 0.03 


14.55 ± 0.03 


15 


.59 


zb 





.03 


15 


52 


± 





.03 


zUUo-U4-ll 


2453837. 


.1 


14. 


.68 ± 0.03 


14.83 ± 0.04 


14.64 ± 0.04 


14.58 ± 0.03 


15 


.53 


zb 





.03 


15 


11 


± 





.04 


2006-04-12 


2453837. 


9 


14 


.95 ± 0.03 


14.95 ± 0.04 


14.68 ± 0.03 


14.49 ± 0.03 


15 


.52 


zb 





.02 


15 


,31 







.03 


2UUo-U4-lo 


24oo8o8. 


.7 


15 


.16 ± 0.03 


15.05 ± 0.03 


14.74 ± 0.03 


14.51 ± 0.02 


15 


.47 


zb 





.02 


15 


,25 


zb 





03 


zUUo-U4-lo 


24ooooy. 





15. 


.27 ± 0.03 


15.17 ± 0.04 


14.80 ± 0.03 


14.52 ± 0.03 


15. 


.43 


i 





.02 


15 


18 


i 

± 





.03 


zUUo-U4-lo 


24ooo4U. 


.8 


15 


.77 ± 0.03 


15.53 ± 0.04 


15.08 ± 0.03 


14.58 ± 0.03 


15 


. 1 1 


i 





.02 


15 


,20 


i 

± 





03 


2006-04-15 


24ooo41. 


.1 


15. 


87 ± 0.04 


15.69 ± 0.05 


15.14 ± 0.04 


14.60 ± 0.03 


15 


.47 


zb 





,03 


15. 


,13 


zb 





03 


o r\ c\ ij n a it 

2006-04-17 


O/l COO /I o 

24ooo42. 


.7 


16. 


,23 ± 0.04 


16.04 ± 0.04 


15.35 ± 0.04 


14.64 ± 0.03 


15 


.42 


zb 





,02 


15. 


,14 


zb 





03 


2006-04-17 


2453842. 


.9 


16 


.30 ± 0.05 


16.22 ± 0.05 


15.40 ± 0.05 


14.66 ± 0.04 


15 


.41 


zb 





.04 


15 


19 


zb 





.05 


ZUUO-U4-1 ( 


-Z40oo4o. 


Q 
O 


10. 


.37 ± 0.04 




15.42 ± 0.04 


14.67 ± 0.03 


lo 


A 1 
.11 


ztz 


n 
u 


-Uo 


it) 




_|_ 

zt 


u 


U4 


2006-04-20 


2453846. 


.1 


17 


01 ± 0.04 


16.91 ± 0.05 


15.91 ± 0.04 


14.94 ± 0.03 


15 


.50 


± 





,03 


15 


,16 


zb 





03 


2006-04-22 


2453847. 


6 


17 


.30 ± 0.05 


17.34 ± 0.07 


16.17 ± 0.04 


15.00 ± 0.03 


15 


.56 


± 





,02 


15 


19 


zb 


0. 


03 


2006-04-24 


2453849. 


.6 


17. 


.73 ± 0.08 


18.13 ± 0.10 


16.65 ± 0.05 


15.18 ± 0.03 


15 


.60 


± 





.03 


15. 


,27 


zb 


0. 


03 


2006-04-26 


2453851. 


.6 


18 


,22 ± 0.09 


>18.8 


17.07 ± 0.07 


15.46 ± 0.03 


15 


.65 


± 





,03 


15. 


23 


zb 


0. 


03 


2006-04-28 


2453853. 


.5 


18 


78 ± 0.10 


>18.8 


17.54 ± 0.07 


15.83 ± 0.03 


15 


.76 


± 





,02 


15. 


23 


zb 





.03 


2006-05-01 


2453857. 







>18.8 


>18.8 


18.29 ± 0.07 


16.31 ± 0.03 


15 


.91 


± 





,02 


15 


22 


zb 


0. 


.02 


2006-05-04 


2453859. 


.7 








18.48 ± 0.11 


16.64 ± 0.04 


15 


.98 


± 





,02 


15 


,22 


zb 





.03 


2006-05-16 


2453871. 


,8 








>18.6 


17.43 ± 0.07 












15. 


.38 


zb 


0. 


02 


2006-05-26 


2453882. 











>18.6 


17.91 ± 0.09 


16 


.64 


± 





,04 


15 


50 


zb 





.03 


2006-05-30 


2453885. 


.9 








>18.6 


>18.2 


16 


.72 


± 





,02 


15 


.56 


± 





02 



comoving frame, subject to the constraints of radia- 
tive and statistical equilibrium (Hillier & Miller 1998). 
The ejecta are assumed to be spherically symmetric and 
chemically homogeneous, a good approximation given 
the very spatially confined regions of continuum and line 
formation (Dessart & Hillier 2005a), and low polariza- 
tion observed in the ejecta of photospheric phase Type 
II SN ejecta (Leonard & Filippenko 2001). The expan- 
sion becomes truly homologous about a week after the 
explosion, but, for simplicity, is assumed homologous at 
all times. We adopt a density distribution characterized 
by a base density po (at the base radius Rq) and a den- 
sity exponent n, with p{R) = p (R /R) n . The value of 
po is chosen so that the Rosseland-mean optical depth 
TRosseiand at i?o is ~50. The hydrogen-rich ejecta are al- 
ways ionized at least once at the base and electron scat- 
tering is the primary source of continuum opacity, while 
the large Hi bound-free cross sections ensure the ther- 
malization of the radiation below the photosphere (see 
§4.2 and Figs 7-8-9 in Dessart & Hillier 2005b). The value 



of n is usually quite large, i.e., ~10, but could be much 
larger at early times, and can be constrained from line 
profile shapes, through the magnitude of the blucshift 
of peak emission of optical P-Cygni profiles (Dessart & 
Hillier 2005a; Blondin et al. 2006; Aldering et al. 2005), 
or through the strength of P-Cygni profiles with respect 
to the continuum. For example, lower values of n lead 
to P-Cygni profiles having stronger and broader absorp- 
tion and emission components, and with a peak emis- 
sion that is closer to line center. The early-time spectra 
of SNe 2005cs and 2006bp, show weak P-Cygni profiles 
with a strongly blucshifted peak emission, and point to- 
wards high values of n. We quote properties of models 
at the photosphere, which corresponds to the location 
where the continuum optical depth, integrated inwards 
from the outer grid radius, is 2/3. 

The luminosity at depth is not known in advance 
unless one uses a consistent hydrodynamical input of 
the ejecta properties, following simultaneously the time- 
dependent evolution of the radiation field, internal en- 
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Fig. 2. — Illustration of the elemental mass fraction distribution 
versus interior mass for the most abundant species in a "reference" 
15 Mq model from Woosley & Heger (2007; a 1.83 Mq neutron star 
is left behind in this explosion). The color coding differentiates 
species: Hydrogen (black), helium (violet), carbon (blue), oxygen 
(turquoise), nitrogen (green), and nickel (red). We also overplot 
the velocity (dotted black line), in units of 10,000km s _1 , to per- 
mit the inference of the elemental stratification in velocity space. 
The surface mass fractions are Xh=0.65, -Xh o =0.34, Xc=0.0014, 
Xn=0.003, -X"o=0.005, and with the adopted mixing, X5e Ni = 1.1 
X 10 — 9 . [See the electronic edition of the Journal for a color version 
of this figure, and see j3]for discussion.] 

ergy, and level populations. Instead, we adjust the base 
(comoving-frame) luminosity so that the synthetic emer- 
gent (observer-frame) luminosity matches approximately 
the observed luminosity at each epoch, for an adopted 
distance and reddening. Usually with our inferred dis- 
tance and reddening, the emergent reddened synthetic 
flux is within a factor of 3 of the observed flux. We 
find that radial/luminosity scaling of that magnitude, 
preserving the ratio L/R 2 produce the same color and 
relative line strengths, line to continuum flux ratios, and 
other essential features of the spectrum, (see Dessart & 
Hillier 2005a for a discussion). 

Reddening is constrained using the Cardelli et al. 
(1988) law for both SNe, assuming R = 3.1. We adopt 
a single value of E(B — V) to obtain satisfactory fits 
to the observed spectra at all epochs. This works best 
when the relative flux calibration is accurate. For SN 
2005cs, we find that the reddening is very small, i.e., 
E(B - V) ~ 0.04 (Baron et al. 2007 infer a value of 
0.035 from a quantitative spectroscopic analysis with 
PHOENIX; Hauschildt & Baron 1999), while for SN 
2006bp, it is quite large, i.e., E(B - V) ~ 0.4 (Ho et 
al. 1997 infer a value of 0.4 to the galaxy host NGC 3953 
from the Ha/H/3 ratio and assuming Case B recombi- 
nation). Swift UVOT observations allow more leverage 
on the estimate of reddening than is permitted by optical 
observations alone by lengthening blueward the observed 
wavelength range. 

To guide our reasoning when interpreting Type II SN 
spectra, in Fig. [2 we illustrate the chemical stratifica- 
tion and velocity distribution versus interior mass for a 
stellar model that is characteristic of the ejecta of Type 
II-P SNe, resulting from the explosion of a red super- 
giant (RSG) progenitor. The figure shows the ejecta of 
a 15 M Q (on the main sequence) progenitor star model 
(ignoring rotation), exploded with a kinetic energy of 
~1 B (=10 51 erg) and with modest (artificial) nickel mix- 



ing (Woosley & Heger 2007) 9 . 

At core collapse, this star has a mass of 
~13M Q (~2M have been lost through a wind during 
the pre-SN evolution), 1.83 M Q is left behind in a 
protoneutron star and the rest escapes to form the SN 
ejecta. Starting from the surface, there is ^7M Q of 
essentially homogeneous material, made up primarily of 
hydrogen (Xh = 0.65) and helium (Xn c = 0.34) - just 
slightly evolved from the primordial, solar metallicity, 
mixture (Grevesse & Sauval 1998). The velocity at the 
base of this hydrogen shell is about 2000 km s _1 . Just 
below is a ~3 M Q shell dominated by helium (Xn c ~ 0.6) 
and oxygen (Xo ~ 0.15), between velocities 1000 to 
2000 km s _1 . Finally, just above the inner mass cut at 
1.83 M Q resides an oxygen-rich shell (Xo ~ 0.6), with 
a ~0.1 mass fraction of helium, carbon, and nickel. 
The radioactive nickel is deep inside, at velocities below 
2000 km s _1 , and separated by ~10 M of material from 
the photosphere at breakout. 

In such (ID) radiation-hydrodynamics simulations, ex- 
plosion parameters are adjusted to obtain the right 
amount of nickel production and mixing, which influ- 
ence the light curve shape or the potential appearance of 
lines from radioactive isotopes at the photosphere (from 
56 Co, for example). For the peculiar Type II SN 1987A, 
this required allowing mixing out to ~4000km s" 1 (Pinto 
& Woosley 1988; Kumagai et al. 1989; Blinnikov et al. 
2000), while for the Type II-P SN 1999em, nickel mix- 
ing out to ^650 km s _1 seems optimal (Utrobin 2007). 
More detailed 2D hydrodynamical simulations of core- 
collapse supernovae explosions, assuming the energy de- 
mands are fulfilled, do not naturally predict nickel at 
large velocities, placing it at ~2500km s~ 1 (Hachisu et 
al. 1990, 1991; Kifonidis et al. 2000, 2003). Observa- 
tionally, there is a range of inferred values and a con- 
servative assumption for Type II-P SNe may be ~1000- 
2000 km s _1 , with 4000 km s _1 representing an extreme 
circumstance. In the above 15 Mq model, 56 Ni is buried 
under about lOM of material, mostly composed of hy- 
drogen and helium, so that any non-thermal excitation 
and ionization of the photosphere by radioactive isotopes 
is delayed until the end of the Plateau phase, when the 
photosphere has receded to these deep and slow-moving 
ejecta layers. In our analysis of the early photosphcric 
phase of Type II-P SN, we neglect any contribution from 
radioactive decay. We expect slow and small variations 
in chemistry during the plateau phase, and, thus, adopt 
chemical homogeneity at a given epoch. This is consis- 
tent with our finding that C and O are only modestly 
enhanced at the photosphere one month after explosion. 
Finally, in such explosions, the density profile follows 
closely a power law distribution with exponent ~10 at 
depth in the hydrogen envelope, but is characterized by 
a much steeper drop-off near the surface layers, that may 
result from the strong radiation pressure in these layers 
at shock breakout (Blinnikov et al. 2000, Shigeyama & 
Nomoto 1990). Such a steepening is not well- modeled 
by Lagrangean simulations, in particular for the shock- 
breakout phase, but our interpretation of our earliest ob- 
servations supports this trend. 

We adopt the surface composition of blue super- 

9 For the CMFGEN spectroscopic modeling that we present here, 
we do not use this ejecta structure at any time. 
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Note. — For each date in our sample of observations, we provide the following CMFGEN model parameters: Base comoving- 
frame luminosity Lcmf,r q and emergent observer-frame luminosity Lobs,j? m (in Lq; note that for the model on the last date, 
the observer-frame luminosity has not fully relaxed to the base condition, a problem that stems from the presence of a steep and 
hard-to-modcl recombination front at the photosphere), photosphcric conditions describing the electron temperature T p hot (in K), 
the radius i? p h t (in 10 14 cm), the velocity Uphot ( m km s — x ), and the density (in 10~ 14 gem -3 ), together with the density exponent 
n characterizing the density law p(r) — p p hot (-Rphot /r) n ■ (See ^3] for a presentation of the modeling approach). a : At is the time 
elapsed since discovery. 



giants analyzed by Crowther ct al. (2006), which is 
in very close agreement with the outer composition 
of the 15 M model described above. By number, we 
adopt H/He = 5, C/He=0.0004, N/He=0.0013, O/He 
= 0.0016, with all metals taken at their solar value 
(Grevesse & Sauval 1998). Given as mass fractions, 
these abundances are: Xn—0.55, Xhc=0.44, Xc=0.0005, 
X N =0.002, Jf o =0.003, and A Fo =0.0013. At later times, 
when in particular 1 7770A gains in strength, we need 
to enhance both the carbon and the oxygen abundances, 
with C/He=0.001, O/He = 0.01, and to reduce the ni- 
trogen abundance, N/He=0.001 to match the data. 

A key asset of CMFGEN is the explicit treatment of 
line blanketing using a detailed description of the atomic 
structure for a large number of species. The iterative 
march to solving many non-LTE level populations to- 
gether with the radiation field comes at great CPU ex- 
pense. We optimize the problem by including only the 
relevant ions of the most important species, although the 
ejecta ionization stratification still requires that at least 
three ionization stages of iron, for example, be included 
to capture the sources of bound-bound and bound-free 
opacity at all depths. Typically, any model crafted for 
very early conditions ( to match the first few, nearly fea- 
tureless spectra of SNe 2005cs and 2006bp, with weak 
Hei5875A) will contain the following species and lev- 
els (the suffixed parenthesis contains the number of full- 
and super-levels; see Hillier & Miller 1998 for details): Hi 
(30,20), He 1(51,40), Hen (30,13), Cn (26,14), Cm (112, 
62), Civ (64,59), Nil (3,3), Nm (26,26), Niv (60,34), 
Nv (67,45), On (3,3), Om (115,79), Oiv (72,53), Ov 
(152,73), Siiv (48,37), Siv (71,33), Fen (309,116), Fern 
(563, 191), Feiv (787,120), Fev (191,47), Fevi (433,44), 
Nim (67,15), and Niiv (242,36). At one week after ex- 
plosion, when He 1 5875A is strong, we take out N iv, O v, 
Si iv, Siv, Fevi, and Niiv, and add Ni (104,54), Oi 
(75,23), Mgn (65,22), Can (77,21), Sin (59, 31), Tin 
(152,37), Tim (206,33), Nin (93,19), and Nim (67,15). 
When hydrogen starts recombining, we put more empha- 
sis on low ionization states and adopt the following model 



atoms: Hi (30,20), Hel(51,40), Ci (63,33), Cn (59,32), 
Ni (104,54), Nai (71,22), Mgn (65,22), Sin (59, 31), 
Oi (75,23), Can (77,21), Fei (136,44), Fen (309,116), 
Fein (477,61), Nin (93,19), and Tin (152,37). Refer- 
ences to the corresponding atomic data sources are given 
in Dessart & Hillier (2005a). 

To summarize, inputs to a CMFGEN model are global 
characteristics of the ejecta such as base (comoving- 
frame) luminosity, radius, chemical composition, velocity 
and density. These parameters are iterated (see Dessart 
& Hillier 2005a for details) until the synthetic spectrum 
matches the observations. Assessment of suitability is 
done by eye. Spectroscopic observations used here are 
accurate in relative flux, while only accurate to within 
10-20% in absolute flux. For spectral fitting, relative 
flux calibration is essential while absolute flux calibra- 
tion is not. When inferring distances, we use the pho- 
tometry to set the absolute flux scale. One important 
result of our study is that the CMFGEN model param- 
eters needed to fit the optical data of SNe 2005cs and 
2006bp during the first 2-3 weeks after explosion also 
naturally explain the Swift UVOT data simultaneously. 
There appears to be a common, photospheric, origin for 
the UV and optical light for both objects: circumstellar 
interaction contributes very little to the observed light in 
these spectral regions (see also Brown et al. 2007). Note 
that SN 2006bp was, however, detected in the X-rays for 
up to 12 days after explosion, suggesting that interaction 
between the ejecta and the circumstellar medium (CSM) 
does occur (Immler et al. 2007), but without any sizable 
effect on the UV and the optical flux. 

The first month after explosion, when the photosphere 
resides in the ionized hydrogen-rich layers of the ejecta, 
is not the best epoch to probe the innermost details of 
the explosion mechanism. However, it is ideal for dis- 
tance determinations based on variants of the Expand- 
ing Photosphere Method: A true continuum exists and 
the effects of line-blanketing are weak in the optical; 
the properties of the ejecta are well determined by the 
model, using both the continuum energy distribution and 
lines from the most abundant ions. The homogeneous 
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Fig. 3. — Comparison between the reddened (E(B-V)=0.04) syn- 
thetic spectrum (red:full; blue: continuum only; Tpi 1C | t = 15750 K) 
and the observations of SN 2005cs on the 30th of June, 2005, 
including CfA observations (black) and Swift photometry (blue 
crosses). We insert a panel that zooms on the optical range to allow 
a better assessment of the fit quality. For completeness, we over- 
lay the Swift UVOT filter passbands in black, (scaled for visibility, 
roughly in proportion to the observed flux), with, from short to 
long wavelengths, uvw2 (A c = 2030A), uvm2 (A c = 223lA), uvwl 
(A c = 2634A), u (A c = 350lA), b (A c = 4329A), v (A c = 5402A), 
where A c is their effective wavelength assuming a Vega-like spec- 
trum (Poole et al. 2007). [See the electro nic e dition of the Journal 
for a color version of this figure, and see £|4.1l for discussion.] 

hydrogen-rich envelope ensures that a well defined pho- 
tosphere exists and that radiation thermalizes at depth, 
allowing the adequate use of the diffusion approxima- 
tion at the model base. There, our non-LTE approach 
naturally predicts LTE conditions, the departure coeffi- 
cients of all levels tending to unity, and at the thermal- 
ization depth, the mean intensity is well represented by a 
blackbody (LTE is, in other words, never enforced). The 
analysis presented here allows the inference of distances 
to Type II-P SNe (see fCDand gO}. 

To be more concise, we present all model parameters in 
Table [6] and [3 but mention in the text only the param- 
eters that are most directly related to the spectroscopic 
evolution, i.e., v p hot, T phot , and n. 

4. SN 2005cs 

SN 2005cs was discovered in the Whirlpool galaxy 
(M51a; NGC5194) by Kloehr et al. (2005) on 28.905 
June 2005 (JD 2453550.4). ft was not detected on 26.89 
June below the limits B > 17.3, V > 17.7, and R >17.6 
(P06), which places a narrow window of two days on 
the explosion date. In this work, we employ observa- 
tions from various sources at 13 epochs, starting on 30 
June 2005, which could be within 3 days of the explo- 
sion. Note that M51a, an SA(s)bc galaxy, also hosted 
the Type Ic SN 19941 (Wheeler et al. 1994; Filippenko 
et al. 1995; Baron et al. 1996). Importantly, Maund 
et al. (2005) and Li et al. (2006) associate the SN 
2005cs event with a 9±l M K0-M4 RSG progenitor star, 
based on the analysis of pre-explosion HST images, com- 
patible with the observed plateau evolution of the light 
curve in the optical bands. SN 2005cs was not detected 
in the X-rays and the upper limit to the X-ray lumi- 
nosity yields a limit to the progenitor mass loss rate of 
M<1 x lO^Moyr" 1 (v w /10)km s" 1 , where v w is the 
asymptotic velocity of the progenitor wind (Brown et 
al. 2007). Model properties are given in Tableland 
corresponding synthetic spectral fits to observations are 
presented in Figs. [3H3 



4.1. The early times: 30th of June until 4th of July 

2005 

We show the first observations (blues crosses: Swift 
UVOT photometry; black curve: optical spectrum of 30 
June 2005 from CfA) of SN 2005cs in Fig. O together 
with synthetic spectra (red: full; blue: continuum only), 
reddened with E(B — V) = 0.04 and scaled by a factor 
1.17. Observations reveal a very hot spectrum, with a 
steep decline of the observed flux longward of ~2300A. 
We identify optical lines of the Balmer series (Ha, H/3, 
and H7), Hel5875A, the Nil multiplets at ~4600A and 
~5400A, and, as suggested by Baron et al. (2007), an 
11 multiplet at ~4600A. Note that the lower ionization 
inferred for the earliest spectra of SN 1999em suggested 
that Nil alone was sufficient to reproduce adequately 
the observed features at both 4600A and 5400A. It is the 
higher ionization in the first spectrum of SN 2005cs that 
calls for this novel 11 identification, more consistent 
with solar CNO chemistry, and slightly different from 
the abundances used in DH06 (see our discussion in $3] 
and Sj7|) . The 4300-5000A region also contains numerous 
He I lines. The observed dip shortward of 2000A stems 
from Fein line-blanketing, rather than dust extinction 
(the continuum-only synthetic spectrum, shown in blue, 
does not show this dip). The individual contributions 
from all important ions on that date are presented in 
Fig. 5, left panel, of Brown et al. (2007). On that date, 
we have Tphot = 15750 K, u p hot = 6880 km s , and n — 
20. Such a steep density distribution is required to bring 
optical lines down to their (weak) observed strength, by 
reducing their emission volumes. This is higher than the 
density exponent adopted in the past, but these are also 
the first analyses of such early-time observations of Type 
II-P SN spectra. A steep density gradient is also used 
to fit the earliest observations of SN 2006bp. We discuss 
this choice and comment on its implications in i j7.ll 

We present P06's observations of SN 2005cs on the 
1st July 2005 in Fig. [U top left panel. The slope of 
the spectrum in the blue appears steeper than predicted 
by the Swift UVOT data points, which are quite closely 
matched by our model (the observed flux in the spec- 
trum, inaccurate below ^4000A, is omitted; Pastorello, 
priv. comm.). The synthetic flux is scaled by a factor 
1.1. The SN ejecta is getting cooler (r pho t = 13350 K) at 
the photosphere, but with a comparable velocity (within 
the uncertainties; v p hot = 6950km s _1 ). The density ex- 
ponent n has been reduced to a value of 12 (the density 
distribution flattens), motivated by the strengthening of 
optical line features, whose identity is the same as that 
on 30 June. The UV flux is weaker, resulting primarily 
from the cooling of the ejecta - the line blanketing is still 
mostly due to Fein. 

Little changes as we move onto the 2nd of July 2005, 
as shown in Fig. [3J top right panel, with a model char- 
acterized by T p hot = 13420 K, v p hot = 6370 km s _1 , 
and n — 12, or a lower ejecta velocity (Dessart & 
Hillier 2005a). This 10% range in expansion velocity 
(^500km s _1 ) may in fact be physical, reflecting the po- 
tential presence of fluid instabilities or turbulence. The 
Swift UVOT photometry would be better fit by a faster 
decline in the UV than our CMFGEN model provides. 
This situation improves on July 3rd, the following day 
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Fig. 4.— Same as Fig. [3] but for the SN 2005cs observations of July 1st (top left; T pho t = 13350 K), 2nd (top right; T phot = 13420 K), 
3rd (bottom left; T phot = 10850 K), and 4th (bottom right; T phot = 9300 K), 2005. A full set of model parameters is given in Table[6] [See 
the electronic edition of the Journal for a color version of this figure, and see ^4. II for discussion.] 



(Fig- SI bottom left panel), where a single model fits both 
the Swift UVOT photometry and the optical spectrum, 
with parameters T p hot = 10850K, v p hot = 6080km s _1 , 
and n = 10. The photosphere is both cooler and slower, 
while the density exponent is now set to 10 at all remain- 
ing epochs. The observations on July 4th resemble those 
of July 3rd in the optical, but the UVOT photometry in- 
dicates a fading of the UV flux. In the bottom right panel 
of Fig. |H we show the model fit to the July 4th observa- 
tions is cooler, with T p hot = 9300 K, u p hot = 5450 km s _1 , 
and n — 10. 

These early observations show a strong UV flux that 
dominates over the optical flux, with low extinction. Op- 
tical lines are Hi, He I, Oil and Nil, with blanketing in 
the UV from Fein. The density distribution is steep, 
flattening over a few days to the canonical n ~ 10 ex- 
ponent predicted for blast wave solutions (Imshennik & 
Nadezhin 1988; Shigeyama & Nomoto 1990, Chevalier 
1982; Ensman k, Burrows 1992, Blinnikov et al. 2000). 
We surmise that the steep density profile may reflect 
work done by the intense UV radiation at and shortly- 
after shock breakout (see H7.ip . Our models support a 
common origin for both the UV and the optical light ob- 
served, with emission from the regions at and above the 
photosphere. No additional contribution to the UV and 
optical fluxes is present from an external source, such as 
an interaction between the SN ejecta and the CSM. 

4.2. The intermediate times: 5th of July until 14th of 
July 2005 

In this epoch, the flux level decreases dramatically in 
the UV and metal lines, mostly of Fen, appear in the 
optical. This is the onset of the hydrogen-recombination 



epoch leading to the "Plateau" phase of such Type II 
SNe. At the same time, Hei5875A weaken, so that the 
feature at ~5900A is due only to NalD at late times. 

These spectral properties are visible in Fig. top 
left panel, in which observations on the 5th of July are 
fitted with a model characterized by Tphot = 8620 K, 
Uphot = 4780km s"\ and n = 10. The Oi7770A line is 
now clearly seen, following the cooling and the recession 
of the photosphere. The SN 2005cs ejecta have lower ve- 
locities than SN 1999em: this allows the clear identifica- 
tion of this line profile, which was severely corrupted by 
the atmospheric extinction in the spectra of SN 1999em 
(DH06). Similarly, while Sin6355A was predicted to 
overlap with Ha in the spectrum of SN 1999em (Figs. 3 
and 5 in Dessart & Hillier 2005a), we clearly observe 
it here as a separate line blueward of the Ha trough 
(see also P06). This line, in combination with uniden- 
tified Nil/On lines at similar Doppler velocities in the 
blue wing of H/3 and Hei5875A, led to speculations on 
the presence of inhomogeneities above the photosphere 
(Leonard et al 2002ab). At early epochs, it seems these 
subtle spectral features can be explained quite naturally 
by numerous lines of On (see above), Nil, Hei, and 
Sin. At later times, the presence of a cool dense shell 
at the outer edge of Type II-P SN ejecta is expected to 
produce real absorption dips in the trough of Ha (and 
Heil0830A; see Chugai et al. 2007). The fit to obser- 
vations on July 6th is shown on the top right panel of 
Fig. [51 with little change from the previous day except for 
a modest decline in temperature (e.g., Tphot — 

8250K). 

By the 9th of July (Fig. middle left panel), the UV 
flux is now an order of magnitude lower than the op- 
tical flux. The CaliH&K doublet appears as a strong 
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Fig. 5.— Same as Fig. El but for the SN 2005cs observations of July 5th (top left; T phot =8620 K), 6th (top right; T phot =8250 K), 9th 
(middle left; T phot =7770 K), 10th (middle right; T phot =7120 K), 11th (bottom left; T phot =7300 K), 12th (bottom right; T phot =6980 K), 
2005 . A full set of model parameters is given in Table [6] [See the electronic edition of the Journal for a color version of this figure, and see 
!]4.2l for discussion.] 
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Fig. 6. — Same as Fig. [3] but this time for P06's observations of SN 2005cs on the 14th of July, 2005, and using a model with enhanced 
(left) or day-one (right) carbon and oxygen abundances. [See the electronic edition of the Journal for a color version of this figure, and see 
q4.2l for discussion.] 
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Fig. 7. — Same as Fig. [3] but this time for the observations of 
SN 2005cs on the 28th of July, 2005. [See the electronic edition 
of the Journal for a color version of this figure, and see £14.31 for 
discussion.] 

P-Cygni line profile at the blue edge of the optical spec- 
trum, strengthening with time, together with the Fe n 
lines present throughout the UV range (as a forest of 
overlapping lines that give an apparently continuous light 
blocking effect - see Fig. 5 in Dessart & Hillier 2005a) and 
up to ~5200A. Mil (shortward of 2500A) and Tin (in 
the 3000A region) also contribute significant line blan- 
keting. The model properties for the 9th of July are 
Tp ho t = 7770K, v phot = 4710km s" 1 , and n = 10. 

Subsequent evolution is quite slow from the 10th 
(Fig. El middle right panel), 11th (Fig. El bottom left 
panel), 12th (Fig. El bottom right panel), and 14th 
(Fig. [6]). Over this 5-day period, the photospheric lay- 
ers cool slowly and modestly (7120 K, 7300 K, 6980 K, 
6700K, respectively), and slow down (4580km a~ l , 
4440 km s"\ 4320 km s" 1 , 3930 km s" 1 ). On the 14th, 
the P06 spectrum extends beyond 9000A and reveals the 
presence of the Can 8500 A multiplet, providing a good 
check that the ionization state of the SN ejecta is well 
modeled. Fitting the Ol7770A line requires a higher 
oxygen abundance (O/Hc=0.01), which we adopt for SN 
2005cs starting on the 11th of July. We show a com- 
parison of the enhanced and standard oxygen abundance 
CMFGEN models in the left and right panels of Fig. El to 
show the effect on the predicted strength of the 1 7770A 
line. Note finally that the Ha line strength (as well as 
H/3, although less discernible) is underestimated by a fac- 
tor of two in both absorption and emission, the more so 
for higher carbon/oxygen abundances (corresponding to 
lower hydrogen mass fractions). 

4.3. 28 July 2005 

The observation of the 28th of July 2005 is the last 
in our SN 2005cs sample. This is about one month af- 
ter the first observation taken and well into the plateau 
phase. We show our fits to these observations in Fig. [7] 
for which the model parameters are T p hot = 6140 K, 
Vphot = 2240 km s _1 , and n — 10. The photosphere 
has now receded to slow-moving ejecta layers close to 
the bottom of the hydrogen shell. 

In practice, given the low luminosity and low velocity 
of this SN envelope, the (outer) hydrogen shell of the 
progenitor may stretch even further down to the slower 
velocity and more C/O enriched layers. Unfortunately, 
the CfA observations do not cover the Oi and Ci lines 
longward of 7500A, which would provide a measure of 



the oxygen and carbon mass fractions. In DH06, and 
by comparison with the surface chemistry of the pro- 
genitor, we inferred significant enhancements in carbon 
and oxygen abundances in the photospheric layers of SN 
1999em at 45 days after explosion. The main line diag- 
nostics are 1 7775A and 1 7984A, while redder transi- 
tions at 8446A and 9264A overlap with Can, Mgn, and 
C I lines, and are therefore less reliable. Useful C I diag- 
nostics are the features at ~9100A, 9400A, and 9600A, 
due to transitions of the (3p-3s) series at 906lA, 9062A, 
9078A, 9088A, 9112A, 9405A, 9603A, 962lA, and 9658A. 
In the near-IR, the Cil.07^m is an additional diagnos- 
tic, although it overlaps with P7, and potentially with 
He 1 1.083 /im. Observations in this wavelength range 
would be especially valuable. 

Observing the C/O abundance change would be im- 
portant and warrants thorough late time spectroscopic 
observations in the photospheric phase of Type II SNe, 
extending at least up to 1/im. A census of lines contribut- 
ing in the optical is shown in Fig. 8 of DH06. Fen blan- 
kets completely whatever UV flux emerges from the ther- 
malization layers below the photosphere. Ti 11 and Fe 11 
leave a complex set of features in the optical, most notice- 
ably resulting in a pronounced dip in the B-band region. 
Na 1 D is reproduced here by enhancing its abundance by 
a factor of four above its solar metallicity value (Grevesse 
& Sauval 1998). However, the predicted line remains 
narrower than observed and time-dependent effects were 
found to lead to a strong and broad line profile without 
any change in sodium abundance. The underestimate of 
the Ha line strength suggests that this is the epoch when 
our steady-state approach departs significantly from the 
more correctly-computed time-dependent conditions in 
the SN ejecta (Dessart & Hillier 2007a,b). 

5. SN 2006bp 

SN 2006bp was discovered in NGC3953, an SB(r)bc 
galaxy, by Koichi Itagaki (Teppo-cho, Yamagata, Japan) 
around 9.60 April 2006 (JD 2453835.1), and located at 
R.A. = Ilh53mn55s.74, Deck = +52°21'09".4 (equinox 
2000.0), which is 62" east and 93" north of the center 
of NGC3953. Quimby et al. (2006) obtained an optical 
spectrum on March 11.11 UT with the 9.2-m Hobby- 
Eberly Telescope (HET) , reporting the observation of a 
blue continuum, a narrow absorption at 592 nm (associ- 
ated with NaiD), and line features at 583 nm (associ- 
ated with He 1 5875 A), and 627 nm (associated with Ha). 
SN 2006bp was not detected in unfiltered ROTSE-IIIb 
images taken on Apr. 9.15 (limiting mag about 16.9), 
but it was found at about 15.3 mag on images taken on 
Apr. 10.15. SN 2006bp was caught promptly after the 
explosion, as confirmed from our analysis of the spec- 
tra. Immler et al. (2007) report the 4.5cr detection of 
X-ray emission at the level of 1.8 ± 0.4 x 10 39 ergs _1 
up until 12 days after discovery, indicating circumstellar 
interaction with a pre-SN RSG wind characterized by 
lO _5 M0yr _1 . This X-ray luminosity is less than 
1% of the UV and optical luminosity in the combined 
UV and optical ranges- and differs from the strongly in- 
teracting Type Iln SNe 10 . Model properties are given 

10 Type Iln spectra are characterized by narrow lines for all ions 
that contribute and most conspicuously Hi and Fell. Examples of 
such Type Iln SNe are 1994W (Chugai et al. 2004), 1995G (Pas- 
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TABLE 7 

Model Characteristics for SN2006bp. 



Julian Date 


Day 


At a 


LcMF,R 


L0BS,fl M ax 


Tphot 


-Rphot 


''phot 


Pphot 

(10- l4 gcm- 3 ) 


n 






Days 


(io 8 L ) 


(K) 


(10 l4 cm) 


(kms" 1 ) 




2453836.6 


2006-04-11 


1.5 


15 


16.3 


20600 


2.386 


14100 


51.7 


50 


2453837.6 


2006-04-12 


2.5 


15 


15.1 


18120 


3.165 


12550 


31.6 


50 


2453842.6 


2006-04-17 


7.5 


9 


6.9 


11800 


5.541 


11300 


8.3 


20 


2453844.6 


2006-04-19 


9.5 


9 


6.6 


11400 


5.868 


11050 


6.5 


16 


2453845.8 


2006-04-20 


10.7 


9 


6.2 


10700 


6.485 


10300 


4.5 


12 


2453846.7 


2006-04-21 


11.6 


9 


6.2 


10700 


6.485 


10300 


4.5 


12 


2453849.6 


2006-04-24 


14.5 


5 


3.6 


9350 


6.345 


7920 


4.2 


10 


2453850.6 


2006-04-25 


15.5 


2.7 


1.8 


8150 


6.580 


8750 


4.0 


10 


2453855.7 


2006-04-30 


20.6 


1.5 


1.1 


7320 


6.000 


6450 


7.3 


10 


2453856.8 


2006-05-01 


21.7 


1.5 


1.1 


7100 


6.147 


6330 


8.8 


10 


2453857.6 


2006-05-02 


22.5 


1.5 


1.1 


7100 


6.147 


6330 


8.8 


10 


2453858.7 


2006-05-03 


23.6 


2.0 


1.4 


6800 


7.26 


5810 


10.1 


10 


2453859.6 


2006-05-04 


24.5 


2.0 


1.4 


6800 


7.26 


5810 


10.1 


10 


2453862.7 


2006-05-07 


27.6 


2.0 


1.7 


6420 


8.90 


4550 


17.9 


10 


2453867.7 


2006-05-12 


32.6 


2.0 


1.7 


6420 


8.90 


4550 


17.9 


10 


2453876.7 


2006-05-21 


41.6 


2.0 


1.8 


6300 


9.28 


4320 


30.3 


10 


2453891.7 


2006-06-05 


56.6 


2.0 


1.8 


6160 


9.73 


3760 


50.7 


10 


2453907.6 


2006-06-21 


72.5 


2.0 


1.9 


6050 


10.18 


3160 


79.4 


10 



Note. — For each date in our sample of observations, we provide the following CMFGEN model parameters: Base comoving- 
frame luminosity Lcmf,h q and emergent observer-frame luminosity LoBS,fl Max (in Lq), photospheric conditions describing the 
electron temperature T p h ot (in K) , the radius R p h.ot ( m 10 14 cm), the velocity u p hot ( m km s - 1 ), and the density (in 10 — 14 g cm -3 ), 
together with the density exponent n characterizing the density law p(r) — p p hot (-Rphot / r ) n ■ (Sec ^3] for a presentation of the 
modeling approach). a : At is the time elapsed since discovery. 




smo flticn flora ioom 

i [A| 

Fig. 8. — Top: Comparison between the reddened (E(B — V) = 
0.4) synthetic flux (red: full; blue: continuum only; Tphot = 
20, 600 K) and SN 2006bp observations on the 11th of April 2006, 
including the optical spectrum of Quimby (2007; black) and Swift 
UVOT fluxes (blue crosses). Note that, despite the larger extinc- 
tion (responsible for the 2200A dip), the UV flux level is compara- 
ble to that observed for SN 2005cs at discovery (shown in Fig- [3J , 
suggesting a much higher ejecta temperature and ionization. Bot- 
tom: Rectified synthetic spectra for the model shown at top, but 
including bound-bound transitions of individual species (see la- 
bel on the right). Apart from the standard Hi and He I lines, 
note the presence of He II 4686A, N III 4001-4099A and N III 4638A, 
C III 2290A and C III 4647A. [See the electroni c ed ition of the Jour- 
nal for a color version of this figure, and see ^15. II for discussion.] 

in Table [7] and corresponding synthetic spectral fits to 

torcllo et al. 2002), 1995N (Fox et al. 2000), 1997eg (Salamanca 
et al. 2002), 1998S (Pooley et al. 2002). Type Iln SNe are strong 
radio and X-ray emitters, associated with the interaction between 
the SN ejecta and the CSM, and with an X-ray luminosity that 
can rival the optical luminosity, as in 1995N (Fox et al. 2000). 



observations are presented in Figs.l8rlT2l 

5.1. First week of evolution of SN 2006bp 

The first spectral observation of SN 2006bp, shown in 
Fig. [5J is impressive in several regards. It has a nearly 
featureless spectrum, with a very steep slope. The dip 
at ~2200A is not caused by line blanketing, as in the 
first spectrum of SN 2005cs (Feiv is the main blocking 
ion in the UV, and it does so blueward of 1800A; its 
effect on the synthetic continuum SED is shown). The 
dip stems from interstellar extinction, with an inferred 
value E(B — V)=0A. We also see, at all epochs, NaiD 
as a narrow absorption feature, indicating significant col- 
umn density in the line of sight to SN 2006bp. Hydro- 
gen B aimer lines are present, but very weak and very 
broad, and limited to Ha, H/3, and H7. Hei5875A is 
just barely identifiable. It is surprising that lines of N 11 
and 11 identified in the first spectra taken of SN 2005cs 
are not seen here. Given the high reddening, fitting the 
steep slope of the SED requires much higher tempera- 
tures and ionization conditions (Fig. I14|). Consequently, 
CMFGEN predicts, although somewhat too strong, the 
presence of Heil4686A (together with a small, ~20% 
contribution from C in 4650A) , a line that is more typ- 
ical of the photosphere/ wind of early- type O and Wolf- 
Rayet stars (Hillier & Miller 1999; Dessart et al. 2000; 
Crowther et al. 2002; Hillier et al. 2003), as well as 
features arising from N in 4001-4099A, Nm4638A, and 
Clll2290A (see bottom panel of Fig. [SJ. In practice, 
we used T p hot = 20,600K, v p hot = 14,100km s _1 , and 
n = 50. Numerous directions in the parameter space 
were explored to reproduce the steepness of the SED, 
the small strength of P-Cygni profiles, and the strong 
blucshift of P-Cygni peak emission. These included vary- 
ing the reddening and temperature, extent of the pho- 
tosphere, ejecta chemistry, and density exponent. The 
most critical was the density exponent. Quimby et al. 
(2007) discuss the presence of weak and narrow lines 
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in this optical spectrum, associated with Heii4200A, 
He II 4686A, and Civ 5808A. They suggest that they em- 
anate from material located just outside the SN ejecta 
that is flash ionized at shock breakout. This gives ad- 
ditional support to our identification here of Hen 4686 A 
emission from the SN ejecta itself. Note the photospheric 
velocity of SN 2005cs inferred from the first observed 
spectrum at 5 days after explosion is only 6880 km s _1 , 
compared to 10 430 km s _1 at a comparable epoch for SN 
1999em; see 2] and DH06). We discuss further the prop- 
erties of SN 2006bp at that age in ffTl 

We could only reproduce this spectrum with CMFGEN 
by increasing the density exponent to 50 (see also 37.1|) . 
Because the density distribution is so steep, line and con- 
tinuum formation coincide and occupy a very confined 
region of space at such early times. As time progresses, 
this region recedes, in the co-moving frame, to deeper 
layers where the density profile is flatter, but the spec- 
trum formation region remains confined, which prevents 
the observer from simultaneously probing the outer re- 
gions. Changes in the exponent that describes the den- 
sity distribution are best inferred through consecutive, 
rather than isolated, spectroscopic observations (Dessart 
& Hillier 2005a). The peak emission of, e.g., the Ha line 
is blueward of its rest wavelength, corresponding to a ve- 
locity shift of ~8000km s _1 initially, but decreases with 
time to reach ~1000km s _1 on the 21st of June. This 
is consistent with a large density exponent (Dessart & 
Hillier 2005a) , although our overestimate of the observed 
Ha-peak blueshift on the 11th of April (Fig. [8]) suggests 
that n may not be as high as 50, and that the weakness 
of optical lines may not be caused entirely by the steep 
density fall-off. 

We show, in the top left panel of Fig. [3J CMFGEN 
fits to observations of SN 2006bp on the 12th of April 
2006. The spectrum looks very similar to that on the 
11th, the fit quality is superior, and the model dif- 
fers in having a lower photospheric temperature and ve- 
locity. The model parameters are T p hot = 18, 120 K, 
Wphot = 12,550 km s _1 , and n = 50. The synthetic UV 
flux also closely matches the flux-equivalent Swift UVOT 
magnitudes, with a nearly featureless spectrum strongly 
extinguished by the high reddening in the line of sight 
of E(B — V) = 0.4. Along a low extinction line of sight, 
the flux at ~2000A would have been about ten times as 
large, reaching a level of ^10~ 13 crgcm~ 2 s _1 A -1 , i.e., 
about a hundred times greater than the optical flux at 
that time! 

The HET observations on the 15th of April 2006 were 
of very low signal and were excluded from our sample. 
For the 17th of April 2006, we show CMFGEN fits in 
the top right panel of Fig. [5J The optical range shows 
hydrogen Balmer lines and He 1 5875A, weak N n features 
in the blue edge of Hei5875A and H/3, as well as a weak 
feature that we attribute to O II lines in the range 4300- 
4600A ( see above discussion for SN 2005cs; Baron et 
al. 2007). The UV flux no longer dominates the opti- 
cal and decreases rapidly from the combined effects of 
cooling (due to radiative losses and expansion) and in- 
creased line-blanketing. The optical lines now appear 
stronger and are well fitted with a density exponent of 
20. The model parameters for the 17th of April 2006 are 
T phot = 11, 800 K, Wp ho t = 11,300 km s"\ and n = 20. 



Two days later, on the 19th of April 2006, the spec- 
trum looks very similar (bottom left panel of Fig. [9]), 
with optical lines somewhat stronger and a diminishing 
UV flux. The model parameters are T p hot = 11, 400 K, 
Wphot = 11, 050 km s" 1 , and n — 16. Note the decreasing 
density exponent from 50 to 16 in just one week, and 
the decreasing photospheric velocity and temperature. 
Because the density distribution flattens, the outer ion- 
ized regions contribute more to the photospheric optical 
depth, shifting it outward to a lower density (Table [7]). 
The same holds on the following day, the 20th of April 
2006 (Fig. [9j bottom right panel) , for which the model 
parameters are Tphot = 10, 700K, w p hot = 10, 300 km s -1 , 
and n = 12. The best-fit density exponent has a value of 
12, and changes little after this time to reach the canon- 
ical value of 10. The spectral evolution is very slow and 
no new lines are identified in the spectrum. N II and O II 
lines have, however, disappeared, following the recombi- 
nation of the ejecta. 

On the 21st of April 2006, we have two sources of ob- 
servations. The data of Quimby et al. (2007) starts 
only at ~4200A but extends out to 10700A, revealing 
a weak 1 7770A and the absence of the Ca II multi- 
plet at 8500A (left panel of Fig. [10]). The Swift UVOT 
data points are well matched by the synthetic UV flux 
distribution, although we overestimate the U-band flux 
somewhat. The model parameters are T p hot = 10, 700 K, 
Wphot = 10, 300km s" 1 , and n = 12. In the right panel of 
Fig. 1 101 we present a similar model with a photospheric 
temperature ~1000K lower, allowing a better match to 
the UV range, but a poorer match to the optical lines. 
This cooler model seems to have the right ionization for 
the UV flux, but has too low an ionization to explain the 
optical. We surmise that this stems from time-dependent 
effects, present at all times, but more visible in spectral 
lines of species (here, helium) that are about to change 
their ionization state in the ejecta (Dessart & Hillier 
2007a, b). This cooler model fits much better the weak 
Fe II 5169A that indicates the start of the recombination 
of Fein to Fen, together with the recombination of hy- 
drogen. So this date closes the early epoch of SN 2006bp 
corresponding to fully ionized conditions in the photo- 
sphere region. 

5.2. The recombination epoch of SN 2006bp 

From this epoch onwards, the ejecta in the photo- 
spheric layers is too cool to radiate a substantial UV flux, 
and line blanketing, due in particular to Fe II lines, blocks 
whatever UV light is produced, re-processing it to longer 
wavelengths. Extinction plays only a secondary role. 

On the 24th of April 2006 (Fig. QTJ top left panel), 
in addition to the strengthening hydrogen Balmer lines 
and the weakening He 1 5875A, we see the appearance of 
spectral features due to Fen lines, around 4400A and 
5200A (the main contributor there is Feii5169A) and 
Mgn4480A. The model parameters are T pho t = 9350 K, 
Wphot = 7920km s -1 , and n = 10. Note that we set the 
density exponent to 10 for this and all following dates. 

On the 25th of April 2006 (Fig. [[TJ top right panel), 
the lines of Fe n strengthen and we see the appearance 
of Si II 6355A, although not as clearly as for SN 2005cs 
due to the higher photospheric velocities. It is not clear 
if the observed dip at around 5800A is due to weakening 
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Fig. 9.— Same as Fig. El but for the SN 2006bp observations of April 12th (top left; T phot = 18120 K), 17th (top right; T pho t = 11800 K), 
19th (bottom left; T phot = 11400 K), and 20th (bottom right; T phot = 10700 K), 2006. A full set of model parameters is given in Table[7] 
[See the electronic edition of the Journal for a color version of this figure, and see ^15. II for discussion.] 
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Fig. 10.— Same as Fig. El but for the SN 2006bp observations on the 21st of April 2006, using the data of Quimby et al. (2007; left) 
or the CfA data (right). The models used in each panel differ mostly in their temperature/ionization, with Tphot = 10700 K (left; Table[7j 
and Tphot = 9700 K (right). [See the electronic edition of the Journal for a color version of this figure, and see ^15. II for discussion.] 



He 1 5875A (as supported by time-dependent calculations 
which predict a higher ejecta ionization sustained over a 
longer time) or due to the strengthening NaiD. In the 
red part of the spectrum, we also see the unambiguous 
presence of 1 7770A and the weak appearance of the 
Caii8500A multiplets, which are not predicted by CM- 
FGEN at this epoch. 

Ha is somewhat overestimated both in emission and 
absorption, which may be related to the assumption of 
steady state or to the constant density exponent adopted 
at a given epoch. The model parameters are T p hot = 
8150K, Up hQ t = 8750km s _1 , and n = 10. 

The spectrum of SN 2006bp on the 30th of April 
2006 (middle left panel of Fig. [TT]) shows evidence of re- 
combination in the ejecta. In particular, Fen lines are 



now clearly visible throughout the F-band, affecting the 
Ha trough (with a small contribution from Si n 6355A) . 
Can lines at 3800A and 8500A are also strong and 
broad. NaiD is the likely contributor to the dip at 
5800A. At this epoch, we employ enhanced oxygen and 
carbon abundance compared to the starting conditions 
(C/He=0.0017 and O/He=0.01 by number), while nitro- 
gen is depleted (N/He=0.0068) and nitrogen lines are not 
seen. In the red part of the spectrum, we observe nu- 
merous Ci lines due to multiplets around 9100A, 9400A, 
and 9630A. We also predict additional contributions from 
O I at 7990A, 8450A, and 9260A, but all much weaker 
than the strong 1 7770A. The model parameters are 
Tp ho t = 7320K, w phot = 6450km s~\ and n = 10. 
For the CfA observations of the 1st and 2nd of May 
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Fig. 11.— Same as Fig. GO but for the SN 2006bp observations of April 24th (top left; T phot =9350 K), 25th (top right; T phot =8150 K), 
30th (middle left; T phot =7320 K), and May 1st (middle right; T phot =7100 K), 2nd (bottom left; T phot =7100 K), and 3rd (bottom right; 
Tphot =6800 K), 2006. A full set of model parameters is given in Table [7] [See the electronic edition of the Journal for a color version of 
this figure, and see £|5.2l for discussion.] 

right), and June 5th (bottom left) and 21st (bottom 
right), 2006. In this elapsed time, the model photo- 
spheric temperature (velocity) decreases steadily from 
6800K (5810km s" 1 ) to 6050K (3160km s" 1 ), while the 
luminosity and the density exponent remain constant (a 
full set of model parameters is given in Table [7]) . 

Throughout this period, the UV flux fades and the 
SED peaks in the optical. The continuum energy dis- 
tribution is very well fitted at all times, becoming more 
and more affected by metal line blanketing, whose ef- 
fect appears first in the optical with Fen lines. Contri- 
butions from Tin lines in the 4200A region strengthen 
with time and cause the B-band magnitude to plum- 
met. By the 21st of June, Fell lines even affect the flux 
in the red part of the spectrum, around 7300A, while 
O I and C I lines continue to contribute numerous narrow 
features in the red. Ca II H&K is overestimated early on, 
while the Can multiplet at 8500A is generally well fit- 
ted. We note the striking evolution of the Ha strength, 



2006bp (middle right and bottom left panels in Fig. [TTj) , 
we employ the same model with Tphot = 7100 K, w p hot = 
6330 km s _1 , and n = 10. The spectral range misses the 
red part of the spectrum where the oxygen, carbon, and 
calcium lines of interest would be found. The spectral 
evolution is slow: we see a slight fading of the flux in 
the blue, with a modest recession of the photosphere, as 
shown by the narrower line profiles. 

We now discuss the late-time observations of SN 
2006bp, an epoch that corresponds, in our steady-state 
modeling with CMFGEN, to conditions where hydrogen 
has recombined at and above the photosphere, and where 
the density of free electrons is low. 

5.3. Late evolution in the photospheric phase of SN 
2006bp 

We present in Fig. rTJ] a montage of synthetic spectra 
fitted to observations of SN 2006bp on May 4th (top 
left), 7th (top right), 12th (middle left), 21st (middle 
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Fig. 12.— Same as Fig. [8] but for the SN 2006bp observations of May 4th (top left; T phot =6800 K), 7th (top right; T phot =6420 K), 
12th (middle left; T phot =6420 K), 21st (middle right; T phot =6300 K), and June 5th (bottom left; T phot =6160 K), and 21st (bottom right; 
Tphot =6050 K), 2006. A full set of model parameters is given in Table [7] [See the electronic edition of the Journal for a color version of 
this figure, and see £|5.2l for discussion.] 



observed with a fixed strength over this 6-week period, 
while steady-state CMFGEN models predict an ever de- 
creasing strength and width (always underestimated) for 
that line (Fig. [T^J). This is caused by time-dependent 
effects, as proposed by Utrobin & Chugai (2005) and 
Dessart & Hillier (2007a, b), and can be remedied by ac- 
counting for the time-dependent terms in the statistical 
and radiative equilibrium equations. Other lines tend to 
be too narrow, in particular NaiD. The good fit to the 
central part of the line profile indicates that the photo- 
spheric velocity is well matched, but the extent of the 
corresponding emission line regions is systematically un- 
derestimated. So, while we provide fits to observations 
at those late epochs, time dependent models would yield 
different parameters and synthetic spectra more compat- 
ible with the observations. Since time dependence does 
not profoundly alter the continuum energy distribution, 
our synthetic fits may still be useful for distance deter- 
minations. It remains to be checked to what extent the 



effects introduced by time dependence affect the correc- 
tion factors used in the EPM, the true magnitude of the 
flux dilution, etc (Dessart & Hillier 2005b). In practice, 
we exclude the observations of SN 2006bp on the 5th and 
on the 21st of June for the determination of its distance. 

6. DISTANCES TO TYPE II-P SNe 

In this section, we compute distances and explosion 
times of SNe 2005cs and 2006bp using our photometric 
observations and synthetic spectra 11 . We use two meth- 
ods, which emerged as the most suitable from the de- 
tailed discussion in Dessart & Hillier (2005b) and DH06. 
The first approach is to apply the Expanding Photo- 
sphere Method (EPM; Kirshner & Kwan 1974; East- 
man & Kirshner 1989; Schmidt et al. 1994; Eastman 

11 Here, we associate the time of explosion with that of shock 
breakout, which may be about a day after core bounce, but corre- 
sponds well to the time when the entire envelope is turned into SN 
ejecta. 
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et al. 1996; Hamuy et al. 2001; Leonard et al. 2002a; 
Elmhamdi et al. 2003), but using correction factors and 
color temperatures calculated directly from our tailored 
models at each epoch (as in DH06). This approach is su- 
perior to using analytical formulae describing the depen- 
dence of correction factors on color temperature, since 
those are computed for models that may not match the 
observed spectrum. This approach avoids the physical 
scatter that stems from such a diversity among SN II. 
The second approach is to use our synthetic spectra to di- 
rectly fit observations, as described in DH06, in a similar 
fashion to the SEAM (Baron et al. 2004). For both, the 
observed magnitudes shown in Table [T] and [2] are used, 
either directly or interpolated linearly between adjacent 
dates to correspond to the date of the spectroscopic ob- 
servation. This does not introduce any sizable error given 
the very slow change (nearly constant from day to day 
within the errors) of J5, V, and / magnitudes in the early 
photospheric phase of Type II-P SNe. 

We have described at length the EPM in Dessart 
& Hillier (2005b), complementing the previous work of 
Eastman et al. (1996), providing some insights into the 
method and its limitations, as well as a new and inde- 
pendently determined set of correction factors based on 
a large database of CMFGEN models. Correction factors 
determined with CMFGEN are comparable to, albeit sys- 
tematically 10-20% larger than, those of Eastman et al., 
although the origin of this small, but important, differ- 
ence remains unclear. In DH06, we applied the EPM to 
SN 1999em and found that, unless the correction factors 
and associated color temperatures are determined with 
tailored models, the EPM distance and the explosion 
time were inconsistent when different bandpass-sets were 
used (amongst BV, BVI, and VI). Since the photo- 
spheric radius is nearly constant across the optical range 
(Dessart & Hillier 2005b) the inferred distance ought to 
be the same whatever optical set is used. 

Another source of uncertainty, although smaller, is 
the inference of the photospheric velocity, w p hot, usu- 
ally approximated by adopting the Doppler velocity at 
maximum absorption in weak and seemingly isolated P- 
Cygni profiles (note that sphericity and projection ef- 
fects, subtleties in the line transfer problem, and di- 
versity in atomic properties lead to a non-trivial rela- 
tion between the photospheric velocity and the velocity 
at maximum absorption in a given line; see Dessart & 
Hillier 2005b for a discussion on this issue). In our ap- 
proach, we adopt the value from our CMFGEN model. 
Although our model fits are satisfactory, they sometimes 
overestimate or underestimate profile widths, matching 
some profiles very well while doing a poorer job on oth- 
ers at the same time. There is likely a genuine range 
of velocities at the photosphere due to fluid instabilities 
and density variations, so that a 10% inaccuracy in the 
photospheric velocity at a given time may in fact exist. 
However, we expect the sense and the magnitude of the 
uncertainty in our inferences to behave randomly. We do 
not consistently overestimate or underestimate v p hot, or 
the flux in a given spectral region, so that observations 
at many times should reduce the uncertainty. 

Reddening is another ingredient requiring special care 
in distance estimates. In our approach, we use synthetic 
SEDs to constrain the reddening E(B — V) within ±0.05. 
The dust properties in other galaxies may not be well re- 



produced by the Cardelli law, which we employ, and this 
introduces an additional uncertainty, which is significant 
for the strongly extinguished SN 2006bp (for such a dis- 
cussion, see, e.g., Elias-Rosa et al. 2006). 

For both objects, we have a good constraint on the ex- 
plosion time from prior optical detection/non-detection. 
However, compact progenitor stars like blue supergiants 
brighten in the optical much more quickly than RSG 
stars with their extended envelopes, whose SED may 
peak in the UV a few days after explosion (the difference 
comes from the energy penalty associated with the ejecta 
expansion). Hence, optical non-detection/detection may 
not constrain the explosion date to better than one day. 
Searching for the UV and X-ray signatures at shock 
breakout would be more valuable in this respect. 

6.1. SN 2005cs 

The high quality spectral fits to the 13 photosphcric- 
phase observations of SN 2005cs presented in the preced- 
ing sections are well suited to determine the distance to 
the SN and its host galaxy M 51a (NGC 5194). 

Focusing first on the EPM, we calculated the correc- 
tion factors and color temperatures for all dates in our 
SN 2005cs sample. As in DH06, we employ three band- 
pass sets, using {B,V}, {B,V,I}, and {V,I}. We log in 
Table [8] the color temperature and associated correction 
factors for each model. Both quantities follow the same 
trend discussed in Dessart & Hillier (2005b) or Eastman 
et al. (1996), with comparable values to those found at 
comparable epochs for SN 1999em (DH06). At late times 
corresponding to the recombination epoch, the correction 
factors rise for {B,V} and {B,V,I}, driven by the strong 
Fell and Till line blanketing in the £>-band. The result- 
ing fading of the B-band magnitude is, thus, not related 
to a variation in the flux-dilution magnitude. Also given 
in Table [8] are the angular sizes resulting from the EPM, 
at each date. Note that these are determined with an 
estimated accuracy of <10%, which is a minimum error 
that accounts for the mismatch between synthetic and 
observed spectra. The SN 2005cs photosphere approxi- 
mately doubles in radius within one month of explosion, a 
comparable increase to that inferred for SN 1999em. Us- 
ing the model photospheric velocity (Table [6]) together 
with these angular sizes, we use the method of least 
squares to fit the equation 9/v p hot = (t — ^explosion )/-£>, 
taking time as the independent variable, and list the re- 
sulting distance and explosion dates in Table [H 

Using either the first 12 dates (to test the effect of a 
shorter time baseline of two weeks, we exclude the last 
date, corresponding to the recombination epoch) or all 
13 dates (encompassing four weeks), we find a distance 
of 8.9±0.5Mpc and an explosion that occurred on JD 
2453547. 6±0. 5, i.e., 2.8 days prior to discovery (~5 days 
prior to the first observation). The agreement between 
the three bandpasses is good. Lengthening the time base- 
line by two weeks by adding the thirteenth date yields 
the same results, which provides further credence to the 
reliability of this result. A plot of these EPM data points 
for the BVI passband, together with a linear fit to the 
data, are shown in Fig. [TS] (top curve). 

The second approach is to fit directly the observations 
with our models. Assuming homologous expansion and 
adopting the model v p \ wt (as in the EPM), we select an 
explosion time to obtain a photospheric radius at each 
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TABLE 8 

EPM QUANTITIES FOR, AND EPM-BASED DISTANCE TO, SN2005CS 



Julian Date 


Day 




Angular Size 




Correction Factor 


Color Tempcrat 


Lire 
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^B V 
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s V I 
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2453552.25 


2005-06-30 


3.006 


3.028 


3.028 


0.562 


0.597 


0.617 


17835 


16813 


16122 


2453553.25 


2005-07-01 


3.778 


3.822 


3.822 


0.496 


0.532 


0.552 


16272 


15280 


14649 


2453554.50 


2005-07-02 


3.645 


3.778 


3.822 


0.496 


0.527 


0.552 


16242 


15400 


14619 
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2005-07-03 
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4.527 


4.527 


0.466 


0.501 


0.527 


14438 


13597 


12935 


2453556.00 


2005-07-04 


5.139 


5.139 


5.104 


0.431 


0.486 


0.526 


13498 


12327 


11456 


2453557.75 


2005-07-05 


5.783 


5.748 


5.692 


0.403 


0.471 


0.518 


12792 


11441 


10525 


2453558.25 


2005-07-06 


5.982 


5.960 


5.872 


0.416 


0.486 


0.542 


12184 


10951 


10020 


2453561.25 


2005-07-09 


6.599 


6.643 


6.709 


0.471 


0.486 


0.491 


10200 


9989 


9899 


2453562.50 


2005-07-10 


6.776 


6.820 


6.864 


0.552 


0.506 


0.491 


8997 


9448 


9689 


2453563.50 


2005-07-11 


6.291 


6.423 


6.555 


0.567 


0.516 


0.486 


9118 


9599 


10080 


2453564.25 


2005-07-12 


6.445 


6.577 


6.754 


0.622 


0.537 


0.491 


8547 


9268 


9929 


2453566.00 


2005-07-14 


6.291 


6.379 


6.423 


0.888 


0.657 


0.542 


6923 


7975 


9178 


2453580.25 


2005-07-28 


7.084 


7.084 


7.216 


1.359 


0.858 


0.642 


5240 


6262 


7525 








BV set 






BVI set 






VI set 





Using First 12 Dates 

D (Mpc) 9.0 ± 0.5 8.9 ± 0.5 8.7 ± 0.6 

texpiosion (JD) 2453547.6 ± 0.5 2453547.6 ± 0.5 2453547.7 ± 0.6 

Using All 13 Dates 

D (Mpc) 9.0 ± 0.5 8.9 ± 0.5 8.7 ± 0.5 

texpiosion (JD) 2453547.6 ± 0.5 2453547.6 ± 0.5 2453547.7 ± 0.5 



Note. — Summary of EPM quantities for SN 2005cs. The quantity 9/v is computed assuming a reddening E(B — V) — 0.04. 
(See text for discussion.) 

obtained with the EPM (here, the quoted error on the 
distance corresponds to the minimum distance-modulus 
scatter between all observations for the chosen time of 
explosion, and is typically 5-10%). We investigate, as 
in DH06 the effect of changing some of the parameters, 
limiting ourselves to the BVI set. Varying the number 
of observations from 13 to 8 moves forward the explosion 
date by about one day while the distance decreases by 
~1.4Mpc. Reducing the time baseline to less than two 
weeks does not give a satisfactory solution: observations 
every three days over the first month is optimal for dis- 
tance determinations based on the SEAM or the EPM. 
Increasing the photospheric velocity at all times by 10% 
(to mimic a systematic overestimate by 10%) places the 
SN further away by 0.9 Mpc, the opposite being true for 
a 10% systematic underestimate of photospheric veloc- 
ity. These uncertainties are gross overestimates since our 
models would not give good fits to the observations at 
all times with such systematic over- or under-estimates. 
Our multi-wavelength multi-epoch modeling to spectro- 
scopic observations of SN 2005cs yields a reddening of 
0.04. A value greater than 0.1 would be impossible to 
accommodate with our models, but values of or 0.1 
would yield acceptable fits, albeit poorer. Keeping the 
same models at all dates, changing the reddening from 
to 0.1 makes the distance vary from 9.4 Mpc down to 
8.2 Mpc. In DH06, we discussed how allowing instead for 
changes in model parameters simultaneously, to adjust to 
the change in reddening, actually reduces the effect of the 
reddening uncertainty. 

Overall, the agreement between the EPM and the 
SEAM-like distances to SN 2005cs give us confidence that 
our inference is accurate (further confidence on such an 
accuracy could come from a Cepheid distance to M 51, 
in analogy to the case of SN 1999em and NGC 1637; see 
Leonard et al. 2003, Baron et al. 2004, DH06). We thus 
adopt a distance of 8.9±0.5Mpc to SN 2005cs and M 51, 
and an explosion date of JD 2453547. 6±0. 5. This in- 



TABLE 9 

Distance to SN2005cs using a direct spectral-fitting 

APPROACH. 
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9 
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8 
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13 
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0.04 
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BV 
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2453548.4° 


8.3 
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Note. — 


Summary o 


f SEAM results 


for SN 2005cs, 


as well 


as tests on 



dependencies with reddening, photospheric velocity, and number of obser- 
vations used in the sample. In this approach, the error on the distance is 
on the order of 5-10% and corresponds to the minimum distance-modulus 
scatter between dates. a : In this case, the distance is obtained by as- 
suming the time of explosion, as shown. JD 2453548.4 corresponds to the 
latest non-detection prior to discovery (P06). (Sec text for discussion.) 

epoch. Each model photospheric radius is then scaled to 
match that corresponding value, together with the syn- 
thetic magnitude. We then calculate the distance mod- 
ulus and the associated standard deviation. We repeat 
this calculation for a range of explosion times, using a fine 
time increment, and select the explosion date that corre- 
sponds to the smallest scatter in distance modulus among 
all dates. The results are given in Table O Using all 13 
dates, a reddening of 0.04, and our model photospheric 
velocity, we find a distance of 8.9±0.7Mpc and an explo- 
sion date of JD 2453547.8, in agreement with the results 
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ferred explosion date is compatible within errors with the 
constraint of 2453549±1 imposed by prior non-detection 
(P06), since even after the explosion has occurred, optical 
brightening is not instant. Adopting a 10% uncertainty 
in the inferred angular size (which corresponds to a min- 
imum error associated with fitting mismatches) leads to 
uncertainties of 0.5 Mpc and 0.5 day for the distance and 
explosion date determined through the EPM. The addi- 
tional uncertainties in photosphcric velocity, reddening, 
and photometry, are small, so that, through our use of 
many observations, our inferred distance should be accu- 
rate to within <10%. 

Our inferred distance is slightly larger than two re- 
cent EPM-like distance measurements to SN 2005cs. Us- 
ing PHOENIX version 14 (Hauschildt & Baron 1999), 
Baron et al. (2007) obtained a SEAM distance to SN 
2005cs of 7.9tol Mpc using all UBVRI magnitudes, but 
at only two epochs. In some sense the level of agree- 
ment is satisfying: PHOENIX and CMFGEN are two 
very distinct codes, and use different assumptions in the 
modeling. Takats & Vinko (2006) used the EPM (with 
analytical correction factors and an estimate of the red- 
dening E(B — V) — 0.1) and the standard candle method 
of Hamuy & Pinto (2002) to determine a distance of 
7.59±1.02Mpc and 6.36±1.3Mpc, respectively. Their 
EPM distance was actually an average of two determi- 
nations — one in which the explosion time was fixed 
(d=8.34±0.31) and one in which the explosion time was 
a free parameter (d=6.84 ± 0.18). The later gave a bet- 
ter fit to the data, but yielded an explosion time 4 days 
later than that inferred directly from the observations. 
Takats & Vinko (2006) find a challenging enhancement 
in the photospheric angular size of about 10-20 over the 
same period, which, in light of our CMFGEN model re- 
sults, would suggest an outward migration, in the La- 
grangean (or co-moving) frame, of the photosphere with 
time. At later times our results, and those of Takats & 
Vinko, appear to be consistent; at early times Takats & 
Vinko estimate much hotter photosphcric temperatures, 
and smaller photospheric radii. 

In general, our inferred distance is larger than most 
other estimates in the literature. Previously, using light- 
curve or spectroscopic modeling of the Type Ic SN 19941, 
Iwamoto et al. (1994) and Baron et al. (1996) obtained 
a distance of 6.9±lMpc and 6±1.9Mpc, respectively, al- 
though the reddening may be uncertain in the line of 
sight to this object. Georgiev et al. (1990) obtained 
a distance of 6.9±0.7Mpc using surface photometry of 
stellar associations. Feldmeier et al. (1997) found a dis- 
tance of 8.4±0.6Mpc using planetary nebulae, revised 
downward by Ciardullo et al. (2002) to 7.6±0.6Mpc 
with an improved reddening of 0.04. Note that this is 
the reddening we infer in the line of sight to SN 2005cs. 
Tonry et al (2001) infer a distance to NGC5195 (M 51b, 
a companion galaxy to M 51) of 7.66±lMpc from sur- 
face brightness fluctuations. Our estimate is compatible, 
within the errors, with that of Tonry et al. (2001), but 
mildly in disagreement with the other studies. Obtaining 
a Cepheid distance to M 51 would help in pinning down 
the distance to that Galaxy and in providing an indepen- 
dent test for the accuracy of our distance estimate. 

6.2. SN 2006bp 



We now apply the methods of the previous section to 
SN 2006bp. To avoid the hydrogen-recombination phase, 
when steady-state CMFGEN models struggle to repro- 
duce the hydrogen Balmer and Na I D lines, we exclude 
the two extreme observations in our sample, those of the 
5th and of the 21st of June 2006, leaving us with 16 ob- 
servations and a time coverage of 4 weeks. 

We present the EPM-inferred quantities in Table [T0l 
giving the angular size, correction factors and color tem- 
peratures for the three sets of bandpasses BV, BVi' , 
and Vi' . Because of the faster expansion of SN 2006bp 
compared to SN 2005cs, the photospheric angular size 
more than triples over the first month of observation. 
At the recombination epoch, the associated front, which 
coincides closely with the photosphere, recedes in the 
co-moving frame at a rate that is comparable to the ex- 
pansion rate of the ejecta. At late times, the angular size 
of the photosphere is, thus, roughly constant. 

Using all 16 observations in our sample and our EPM- 
based technique leads to a distance of 17.5±0.6Mpc for 
the BVi' set, only slightly varying with bandpass set 
(17.9±0.7Mpc for the BV set and 17.4±0.6Mpc for the 
Vi' set; Table ITUf . Including only the first 13 dates (to 
assess the effect of a shorter time baseline of one month) 
leads to a higher distance by 0.2-0.3 Mpc in each band 
and a larger error of ±0.8 Mpc. The explosion date is also 
in good agreement for either 13 or 16 dates in our sample, 
with a Julian date of 2453833. 4±0.4day. A plot of these 
EPM data points for the BVi' passband, together with 
a linear fit to the data, are shown in Fig. [13] (bottom 
curve) . 

We also use our synthetic spectra and magnitudes to 
constrain directly the distance, in a fashion similar to the 
SEAM. We find a distance in agreement with the EPM 
inference, with 17.1 Mpc for the BVi' set (Table [TT|). 
Uncertainties are larger than the scatter of values across 
bandpasses, but should remain at the 5-10% level, in 
particular due to the large number of observations. Er- 
rors are unlikely to be purely statistical, but they have 
some level of randomness as the deviations from a per- 
fect match of synthetic to observed spectra are very dif- 
ferent for individual epochs (blue featureless spectrum 
at early time to strongly blanketed spectrum at the re- 
combination epoch), different spectral regions (the blue 
and red wavelength regions reflect different contributions 
from different elements that we model with variable lev- 
els of success). Reducing the number of observations to 
9 dates increases the distance estimate by 1.7 Mpc. Un- 
certainties in reddening and photospheric velocity have 
a similar effect to that described for SN 2005cs (JI]) and 
SN 1999em (DH06). 

To conclude, we estimate a distance to SN 2006bp and 
its galaxy host NGC 3953 of 17.5±0.8Mpc together with 
an explosion date of 2453833. 4±0.4day only 1.7 days af- 
ter discovery (~3 days prior to the first observation). 
This makes it the earliest-discovered Type II-P SN and 
gives further support to the novel properties of the SN 
ejecta inferred from the first observed spectra (steep den- 
sity distribution, highly ionized ejecta with Hen4686A 
line emission). This distance estimate is in close match 
to the value of 17 Mpc obtained by Tully (1988). This 
explosion date is also in agreement with the date of April 
9th, inferred by Quimby et al. (2007) from fits to optical 
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TABLE 10 

EPM QUANTITIES FOR, AND EPM-BASED DISTANCE TO, SN2006BP 
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Using First 13 Dates 
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(Mpc) 








18.1±0.8Mpc 






17.8±0.8Mpc 






17.7±0.8Mpc 




t explosion ( JD) 








2453833.2±0.4 






2453833. 3±0.4 






2453833. 3±0.4 





Using All 16 Dates 



D (Mpc) 17.9±0.7Mpc 17.5±0.6Mpc 17.4±0.6Mpc 

toxpiosion (JD) 2453833.3±0.3 2453833. 4±0.3 2453833. 4±0.3 



Note. — Summary of EPM quantities for SN 2006bp. The quantity 9/v is computed assuming a reddening E(B — V) — 0.4. 
(See text for discussion.) 



TABLE 11 

Distance to SN2006bp using a direct spectral-fitting 
approach. 
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Note. — Summary of SEAM results for SN 2006bp, as well as tests on 
dependencies with reddening, photospheric velocity, and number of obser- 
vations used in the sample. In this approach, the error on the distance is 
of the order of 5-10%, and corresponds to the minimum distance-modulus 
scatter between dates. a : In this case, the distance is obtained by assum- 
ing the time of explosion, as shown. (See text for discussion.) 

ROTSE-III light curve (note however that the SN may 
remain dim in the optical for many hours after breakout). 

7. DISCUSSION 

7.1. Ejecta properties a few days after shock breakout 

We have presented quantitative spectroscopic anal- 
yses of two Type II-P SNe, using observations trig- 
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Fig. 13. — Linear fit (dotted line) to the BVI data shown 
in Table |8| for SN 2005cs (using 13 dates) and to the BVi' data 
shown in Table [Tul for SN2006bp (using 16 dates). The reciprocal 
of the slope is the distance D, i.e. 6/u p h t = (* — ^explosion)/^ 
where 9 = R p h t/D. The origin for the time is the inferred time of 
explosion (more specifically the shock breakout time) . The vertical 
bars indicate a ±10% error. 

gered at an exceptionally early time after explosion 
(3-5 days), providing a glimpse into the ejecta struc- 
ture relatively soon after shock breakout. We infer a 
very high-temperature/very high ionization in the corre- 
sponding photospheric layers and in particular identify 
the presence of broad, ejecta-related, He n 4686A in the 
first two observations of SN 2006bp (narrow emission in 
He II 4686 A is also identified by Quimby et al. (2007), 
related to material outside of the SN ejecta, which also 
testifies for the high temperature associated with shock 
breakout). In Fig. [TJ] we show the ionization level for 
H, He, C, N, O, and Fe in the ejecta computed with 
CMFGEN and for models that fit the first spectroscopic 
and photometric observations of SN 2006bp (left panel) , 
SN 2005cs (middle panel), and SN 1999em (right panel). 
These correspond, in the same sequence, to 3, 5, and 6 
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days after explosion, as determined in this work and in 
DH06 12 . One clearly sees in this figure the evolution to 
lower ionization through this sequence, following directly 
the elapsed time since explosion. Oxygen changes from 
being twice- to singly-ionized at mid- heights, from SN 
2006bp to SN 1999em. Helium is singly ionized through- 
out the SN 1999em ejecta, but twice ionized in the in- 
ner and outer regions of the SN 2006bp ejecta. This 
trend towards much higher ionization at earlier times 
supports the identification of He II 4686A emission in the 
SN 2006bp ejecta, as well as O n lines for both SN 2006bp 
(at about 5 days after explosion) and SN 2005cs (in the 
discovery spectrum). On lines were not detected un- 
ambiguously in the first spectrum taken of SN 1999em 
because it was taken too late, about 1-2 days later than 
SN 2005cs relative to the explosion time. 

Moreover, at such early times, adopting the nitrogen 
abundance of DH06 now overestimates the strength of 
N II lines and underestimates the strength of O II lines, 
justifying the adoption in this work of CNO abundances 
closer to solar, i.e. C/He=0.0004, N/He=0.0013, O/He 
= 0.0016, in agreement with the surface chemistry of 
blue supergiants analyzed by Crowther et al. (2006). By 
comparison, DH06 used C/He=0.00017, N/He=0.0068, 
O/He = 0.0001 (all given by number; for both works we 
use H/Hc=5). We thus concur with the findings of Baron 
et al. (2007) on the necessary CNO chemistry to repro- 
duce the observations of O II and N II lines in the early 
spectra of SN 2005cs (as well as in SN 2006bp), with a 
depletion of carbon and oxygen compared to solar that 
is modest rather than severe. In general, an accurate 
determination of CNO abundances in the fast-expanding 
photospheres of Type II-P SNe is difficult since the cor- 
responding lines are weak and overlap with lines of other 
species such as He I. Having recourse to quantitative 
spectroscopy of the emergent light from their supergiant 
progenitors, which can be found in number and nearby, 
seems a more viable route to constraining the surface 
chemistry of such pre-SN massive stars. 

In both objects, the near featureless spectral appear- 
ance over the first few days requires the choice of a very 
high power-law density exponent of 20-50 to dwarf the 
emission volume of all optical lines, as well as very high 
ejecta velocities to reproduce the large widths of emission 
line features. This is admittedly well above the value 
n = 10 for the propagation of a shock wave through 
the envelope of a massive star progenitor, both obtained 
theoretically (Imshennik & Nadezhin 1988; Shigeyama 
& Nomoto 1990, Chevalier 1982; Ensman &; Burrows 
1992, Blinnikov et al. 2000) and inferred observation- 
ally (see, e.g., Eastman & Kirshner 1989; Schmutz et al. 
1990; Baron et al. 2000; Dessart & Hillier 2005a,2006a). 
Nonetheless, for the outermost layers of the progenitor 
envelope, a much steeper density profile is expected, with 
values in the vicinity of a few tens (Blinnikov et al. 2000, 
and in particular their Fig. 1), and, under certain condi- 
tions, a small amount of mass at the progenitor surface 
may even reach relativistic speeds as the shock breaks 
out (Matzner & McKee 1999). 

12 The maximum radius differs between simulations, resulting 
from the different density exponents, but in all cases is such that 
the radiation is freely-streaming at the outer grid radius (optically- 
thin conditions) 



Unfortunately, the radiation-hydrodynamics codes 
used to compute the ejecta structures of core-collapse 
SN are Lagrangean, and are thus not well designed to 
describe the dynamics and the structure of the low den- 
sity surface layers. In those fully-ionized regions, the 
luminosity is super-Eddington by up to a few orders of 
magnitude, even when accounting only for the electron- 
scattering opacity. When accounting for the additional 
contribution of optically thick and thin lines, we find 
that the radiation pressure can be very substantial, fol- 
lowing the combination of a huge luminosity, a low den- 
sity, and a large velocity gradient (Castor, Abbott, & 
Klein 1975). In the CMFGEN model of SN 2006bp on 30 
June 2006, the radiative acceleration in the layers above 
the photosphere reaches a few thousand cms -2 , which 
combined with a velocity of a few 10000 km s -1 and a 
length scale of 10 14 cm lead to a velocity gain of few 
tens of percent (assuming the acceleration remains con- 
stant). Sooner after shock breakout, the radiative accel- 
eration could be even larger. Overall, we expect that the 
outer optically-thin layers would adopt a steeper den- 
sity distribution and a faster-than-linear velocity distri- 
bution with radius. We plan to investigate these issues 
more quantitatively in the future, but at present, it seems 
that the earliest observations of SNe 2005cs and 2006bp 
support the very steep density distribution obtained in 
radiation-hydrodynamics simulations of core-collapse SN 
ejecta soon after shock breakout. 

Before closing this section, we note that SN 1999gi is 
another Type II-P SN discovered very early, perhaps only 
~4+3° days after explosion (Leonard et al. 2002b). In- 
terestingly, Leonard et al. report the presence of a full 
P-Cygni profile shifted from the rest wavelength of H/3 
by —30000km s _1 in the first spectrum of SN 1999gi, and 
speculate on the potential association with ejecta inho- 
mogeneities. In the context of a density inhomogeneity 
far above the photosphere, a blueshifted absorption is 
a possibility, but a P-Cygni profile entirely blueshifted, 
and by that amount, is not (see, however, Dessart & 
Hillier 2005a). We reproduce their observation in Fig.[T5l 
and overplot the first spectrum of SN 2006bp (red line). 
Given the close correspondence between the two spectra 
and the likely early detection of SN 1999gi, we propose 
that, in fact, the broad line feature just blueward of H/3 
in the spectrum of SN 1999gi is due to Hen 4686A, as for 
SN 2006bp. 

7.2. Energetics 

Having inferred the reddening and the distance to SNe 
2005cs and 2006bp, we can deduce the absolute mag- 
nitude for each event at all epochs studied. Together 
with results presented for SN 1999em in DH06, we show 
in the top-left panel of Fig. [T§] the evolution with time 
since explosion of the absolute magnitude in B (black), V 
(blue), and / (red), for SN 2005cs (solid line), SN 2006bp 
(dashed line; i 1 is shown instead of /), and SN 1999em 
(dotted line). Note that the time of explosion is com- 
puted together with the distance, rather than estimated 
based on prior non-detection of the SN on pre-explosion 
images of the corresponding host-galaxy 13 . While the 

13 Our inferences are, however, compatible with these alternate 
estimates, in particular for SN 2006bp for which the explosion time 
is very tightly constrained by a non-detection on April 9.15 and 
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Fig. 14. — Radial variation of the ionization state of H (black), He (blue), C (turquoise), N (green), O (yellow), and Fe (red) in the ejecta 
of SN 2006bp (left; ij =2.2xl0 14 cm), SN 2005cs (middle; ij =1.6xl0 14 cm), and SN 1999em (right; R =4.176xl0 14 cm) based on the 
modeling of the first spectroscopic observation (which corresponds to 3 ([[5]and Fig. [8}, 5 ([|4]and Fig.0, and 6 (DH06 and their Fig. 1) 
days after explosion) and limited to the regions exterior to the layer where the inward integrated Rosseland optical depth is ~50 on that 
day. Note the systematic decrease of the ejecta ionization state as we go from SN 2006bp, to SN 2005cs, and to SN 1999em, following 
the trend of elapsed time since explosion. This clarifies the identified emission in Hen4686A line in the first spectrum of SN 2006bp. 
Note also the rise in ionization at large distances (together with a rise of the temperature, not shown here), a non- LTE effect in these 
scattering-dominated environments. [See the electronic edition of the Journal for a color version of this figure, and see £|7.1l for discussion.] 
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Fig. 15. — Comparison between the first spectroscopic observa- 
tion of SN 2006bp (11th of April 2006) and of SN 1999gi (10 De- 
cember 1999; Leonard et al. 2002b) showing the observed flux ver- 
sus the Doppler velocity, with respect to the rest wavelength of H/3 
(Note that the inferred reddening of SN 1999gi is E(B - V) = 0.21, 
so about half what we infer for SN 2006bp; Leonard et al. 2002b). 
The P-Cygni profile seen in the SN 2006bp spectrum is identified as 
Hell4686A (Fig. [Sj , and, given the similarity in profile shape and 
position, suggests a similar identification for SN 1999gi. [See the 
electronic edition of the Journal for a color version of this figure, 
and see ^7. II for discussion.] 

spectroscopic and photometric evolution of three events 
are comparable in all spectral regions at a given time 
after explosion, they display a sizable range of intrinsic 
brightness. SN 2005cs is about two magnitudes fainter 
than SN 2006bp, about a factor of six in (bolometric) lu- 
minosity, with SN 1999em somewhat less luminous than 
SN 2006bp. SN 2006bp was detected earliest of all three 
and is in fact the first Type II-P to be caught so soon 
after explosion, with a ^2-day delay at most (note that 
1987A and 1993 J are not Type II-P SNe). 

In the top-right panel of Fig. Ql)l we show the corre- 
sponding evolution of the model photosphcric velocity, 
for SNe 2005cs, 1999em, and 2006bp. The above trend is 
maintained, i.e., larger intrinsic brightness is associated 
with larger ejecta velocity. Note that if the bolometric 
luminosity scaled with the square of the ejecta velocity, 

detection on April 10.15. 



a brightening by 1.5 mag would result from a factor of 2 
enhancement in velocity, which is compatible with our in- 
ferences for SNe 2005cs and 2006bp (top row of Fig.HBJ). 
Since we adopt similar ejecta structures (after one week 
past explosion), this suggests that the bolometric lumi- 
nosity correlates with the ejecta kinetic energy in the 
hydrogen shell, as pointed out in the past by, for exam- 
ple, Falk & Arnett (1977), Litvinova & Nadezhin (1983), 
Hamuy (2003), Zampieri et al. (2003), Pastorello et al. 
(2004). The late time brightness is not controlled by 
the energy deposited at shock passage, but rather stems 
from the amount and location of the unstable nickel at 
the start of homologous expansion, so the above trend 
may not carry over until much later times, although it 
might, but for different reasons (higher energy explosions 
may systematically produce more 56 Ni). 

Similarly, and for completeness, we include in the 
bottom panels of Fig. [16] the temporal evolution of 
the photospheric radius (left) and temperature (right). 
Note the smaller extent of the slowly-expanding ejecta 
of SN2005cs, the converging photospheric temperatures, 
and the near-constant photospheric radii at the recom- 
bination epoch for all three SNe, tightly linked to the 
Plateau appearance of the visual light-curves at such 
times. 

7.3. A common origin for the UV and the optical flux? 

Our modeling approach has been first to fit the opti- 
cal range, using both the overall shape of the SED and 
spectral lines formed in the SN ejecta. As a consistency 
check, as well as an additional constraint on the result- 
ing model parameters and the choice of reddening, we 
compared the synthetic spectrum to the observed Swift 
UVOT values. For both SN 2005cs and SN 2006bp, we 
obtain very satisfactory fits and we thus conclude that 
the observed UV and optical light emerges from the pho- 
tosphcric layers of the SN ejecta, over the first 1-2 months 
after explosion (starting 3-5 days afterwards). This re- 
sult is unambiguous, and, thus, no sizable contribution 
to the UV and optical fluxes stems from an interaction 
with the CSM. 



24 



•a 
3 

'S 
m 

CO 

s 



-17 : 



18 r 



-15 - 



-14 




Dashed: 06bp 
Dotted: 99em 
Solid: 05cs 
BVI i 



10 20 30 40 

Time after Explosion [days] 



50 





14000 




12000 


l 




1 


10000 


H 




S 


8000 






4 


6000 


> 


4000 




2000 









i 




Dashed: 06bp 
Dotted: 99em 
Solid: 05cs 



10 20 30 40 

Time after Explosion [days] 



50 




Fig. 16.— Top Left: Evolution of the absolute magnitude in the B (black), V (blue), and / (red) bands for SN 2005cs (solid), SN 2006bp 
(dashed; i' is shown instead of I), and SN 1999em (dotted), with respect to the inferred time of explosion (JD 2453547.6 for SN 2005cs; 
JD 2453833.4 for SN 2006bp; JD 2451474.8 for SN 1999em), and the inferred distance of 8.9 Mpc and reddening of 0.04 for SN 2005cs, 
11.5 Mpc and 0.1 for SN 1999cm, and 17.5 Mpc and 0.4 for SN 2006bp. Note how the SN 2005cs is intrinsically underluminous compared 
to SN 1999em, by about 1.5 mag equivalent to a factor of four in luminosity (a valid comparison since their SED is comparable at a given 
time after explosion). SN 2006bp is somewhat more luminous than SN 1999em, by about 0.5 mag in each bandpass. Crosses refer to the 
times of the observations that we use here, and given with respect to the inferred time of explosion obtained in this work and in DH06 
for SN 1999em. Top Right: Same as left, but this time for the photospheric velocity. Notice how the relationship in absolute magnitude 
carries over for the photospheric velocity (uncertain at the 10% level), at the same time after explosion (note that the three ejecta have 
a similar density structure beyond one week after explosion). The larger the kinetic energy, the larger the bolometric luminosity of the 
ejecta. Bottom: Same as the top right panel, but for the photospheric radius (left) and temperature (right). [See the electronic edition of 
the Journal for a color version of this figure, and see >]7.2l for discussion.] 



A further evidence for the photospheric origin of the 
UV flux is the identification of the Mgn2800A line as 
a broad P-Cygni profile in Swift UVOT grism spectra 
(Brown et al. 2007), with a width compatible with the 
photospheric velocity at the corresponding epoch. So, 
while circumstellar interaction may be important in some 
Type II SNe (the Iln events), it seems to contribute neg- 
ligibly to the UV and optical light we received from SN 
2005cs and SN 2006bp, and more generally in agreement 
with the weak radio and X-ray emission from Type II- P 
SNe (see, e.g., Pooley et al. 2002; Chevalier et al. 2006). 

7.4. Comparison with previous work 

The recent spectroscopic analysis by Baron et al. 
(2007) of SN 2005cs at two epochs in the photospheric 
phase agrees in many respects with the findings of the 
present study for that object and, thus, gives some cre- 
dence to the Type II SN spectroscopic analyzes per- 
formed with PHOENIX (Hauschildt & Baron 1999) and 
CMFGEN (Hillier & Miller 1998; Dessart & Hillier 
2005a). This is particularly satisfying since PHOENIX 
and CMFGEN are two very distinct codes, and use differ- 
ent assumptions in the modeling. A potential important 
difference in assumptions is our neglect of energy con- 
tribution from radioactive decay throughout this work. 
By contrast with Baron et al. (2007), we find no strong 



evidence for nickel mixing and non-thermal energy depo- 
sition at the photosphere over the first month after ex- 
plosion (in this work, we obtain good fits to observations 
prior to the hydrogen-recombination epoch without in- 
voking this additional source of energy). The ejecta ion- 
ization structure computed by the two codes may thus 
differ. 

In this work, we infer a reddening E(B — V) ~ 0.04, in 
agreement with Baron et al. (2007) who obtain a value 
of 0.035. We also concur with the findings of Baron 
et al. (2007) on the necessary CNO chemistry to re- 
produce the observations of Oil and Nil lines in the 
early spectra of SN 2005cs (as well as in SN 2006bp), 
with a depletion of carbon and oxygen compared to so- 
lar that is modest rather than severe. As discussed in 
^7.11 this constraint can be placed for SN 2005cs, by 
contrast with, e.g., SN 1999em, mostly because it was 
observed at an earlier time when the ionization con- 
ditions at the photosphere are ripe for On line emis- 
sion. The identification of On and Nil lines can ex- 
plain most, if not all, the mysterious kinks seen in the 
blue wing of H/3 or Hei5875A (see, e.g., Leonard et al. 
2002a, b who reported such kinks in the spectra of SNe 
1999em and 1999gi). Here, the abundances adopted at 
early times for both SNe 2005cs and 2006bp are, given 
as mass fractions, Ah=0.55, Ah g =0.44, Ac=0.0005, 
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A N =0.002, A o =0.003, and A Fe =0.0013. At later times, 
when in particular 1 7770A gains in strength, we en- 
hance both the carbon and the oxygen abundances, with 
C/He=0.001, O/He = 0.01, and reduce the nitrogen 
abundance, N/Hc=0.001. Again, observations support 
the same ejecta chemistry at the photosphere at cor- 
responding epochs after explosion for both SNe. We 
are however reserved concerning the sodium abundance, 
which we predicted, based on steady-state CMFGEN 
models, to be enhanced by a factor of four over its pri- 
mordial value in SN 1999em (DH06). Comparable en- 
hancements would be needed here for 2005cs and 2006bp, 
and, indeed, Baron et al. (2007) predict a sodium en- 
hancement of a factor of 10. Time- dependent CMFGEN 
models performed with a primordial sodium abundance, 
however, suggest that Nai D may remain as strong and 
broad as observed (Dessart & Hillier 2007a, b); further 
work is needed to quantify the effect of time-dependence 
on abundance determinations. 

For the observations of the 14th of July 2005, Baron 
et al. (2007) find at the photosphere a temperature of 
6000 K, a velocity of 4000 km s -1 and a density exponent 
of 8. In the same order and for that date, we find 6700 K, 
3930 kms" 1 , and 10, in agreement with their values. For 
the observations of the 31st of July 2005, Baron et al. 
find at the photosphere a temperature of 5500 K, a ve- 
locity of 6000kms _1 and a density exponent of 12 (see 
their §6). In this work, for the observations of the 28th of 
July 2005, just three days before (note that observations 
do not change visibly over just a few days at such late 
times), we find in the same order 6140 K, 2240kms _1 and 
a density exponent of 10. While we fit the CfA obser- 
vations reasonably well with this parameter set on that 
date (Fig[7j), Baron et al. require a photospheric velocity 
nearly three times as large to do so, and also larger than 
the value on Day 17 (in fact nearly as large as on the first 
spectrum taken, for which we use 6880 km s -1 ). Upon in- 
spection of their Figs. 12-13, it appears that Ha peaks 
at ~6480A, while the observations they use, as originally 
published in P06, and our (restframe) CfA observations 
show a peak at 6550A. This difference of 70A corresponds 
to a velocity shift of 3200 km s -1 at this wavelength and 
suggests that their extravagant proposition of an out- 
ward moving photosphere between Day 17 and Day 34 
(see their §6) may stem solely from a problem with the 
wavelength of their flux datapoints. 

Baron et al. (2007) obtained a SEAM distance to SN 
2005cs of 7.9±a 7 6 Mpc using all UBVRI magnitudes, but 
at only two epochs, which is in agreement with our esti- 
mate of 8.9±0.5Mpc using BVI and 12 or 13 epochs. 

Although the observed spectra of SNe 2005cs and 
2006bp presented in this work are different in many ways, 
these reflect mostly the difference in the ejecta expansion 
rate (small for 2005cs), the reddening in the line of sight 
(large for 2006bp), and the elapsed time between explo- 
sion and discovery (very short for 2006bp). As shown 
in Fig. [TBI both SNe evolve in a similar fashion, 2005cs 
being merely undcrluminous and boasting a lower ejecta 
kinetic energy. 

8. CONCLUSION 

We have presented a multi-epoch multi-wavelength 
quantitative spectroscopic analysis of the two Type II-P 
SN events 2005cs and 2006bp, over the first two months 



after explosion, and have reached similar conclusions on 
the ejecta properties to those for SN 1999em (DH06). 
The progenitor stars have a composition compatible with 
a BSG/RSG progenitor. The spectral evolution shows a 
reddening of the SED with time, with a flux peaking in 
the UV at very early times, but progressively shifting to 
longer wavelengths. The photosphere recedes to deeper 
layers with time, reaching mass shells moving at about 
2000-3000 km s _1 about 6 weeks after explosion and hav- 
ing enhanced carbon and oxygen, and reduced nitrogen, 
abundances, compared to the surface layers. We find 
a density distribution that is very steep at early times, 
and we surmise that this could be an imprint of the shock 
breakout epoch. From about a week until six weeks after 
explosion, we find a more standard and constant density 
exponent of ten for the corresponding power law distri- 
bution, compared to a value of 50 (for SN 2006bp) and 
20 (for SN 2005cs) used to model the first observation. 
Swift UVOT photometry and ground-based optical spec- 
troscopy and photometry support a common origin for 
the UV and the optical light, so that we anticipate that 
if CSM interaction does occur, it does not contribute sig- 
nificantly to the observed UV and optical fluxes. 

Using our analysis, we infer the distance and the ex- 
plosion date of each SN. We find, in that order, JD 
2453547. 6±0. 5 days (2.8 days prior to discovery) and 
8.9±0.5Mpc for SN 2005cs, and JD 2453833.4±0.4days 
(1.7 days prior to discovery) and 17.5±0.8Mpc for SN 
2006bp. The large number of observations, the quality 
of our fits, the agreement between the two methods fol- 
lowed suggests a very high level of internal consistency 
of the distances determined here with Type II SNe. The 
actual level of accuracy on such determined distances, 
which we estimate to be on the order of 10% (related 
to the spectral fitting accuracy), would be best assessed 
by confrontation with alternate methods (Cepheids for 
example, although all alternate methods have their own 
inaccuracies, claimed and otherwise), the determination 
of distances to a galaxy host with well-recorded multi- 
ple Type II SN events (the typical Type II SN rate of 2 
per century per galaxy makes this proposition somewhat 
unrealistic), a Galactic SN that could also be resolved 
with an interferometer or a large- aperture telescope. Our 
perspective on the uncertainty of our distance measure- 
ments is that, were we to increase the number of Type 
II SNe suitable for analysis in M 51 to 10 or 100, the 
scatter in the distance determined with our technique 
for that sample would not be larger than 10%. But this 
would still be no proof that the resulting average dis- 
tance is indeed the right one since enlarging the sample 
would only allow for an internal consistency check on 
our method. In the short term, a Cepheid distance to 
the galaxy host of 2005cs and 2006bp is highly desirable, 
and would provide an independent check on the distance 
computed with Type II-P SNe, as was done with success 
in the past with SN 1999em (Leonard et al. 2003; Baron 
et al. 2004; DH06), giving further support to the poten- 
tial use of such objects for distance determinations in the 
Universe. 

This work fits within an ambitious project to build a 
general understanding of SN ejecta, to infer properties of 
the progenitor star (with ramifications for massive star 
evolution) and the explosion mechanisms at the origin of 
the event. A major ongoing theoretical development is 
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the treatment of time-dependent effects both in the ra- 
diation field and in the level populations to allow for a 
direct modeling of the evolution over months after shock 
breakout, through the photospheric and into the nebular 
phases, and based on a range of hydrodynamical models 
of the explosion. In parallel, high-quality multi-epoch 
multi-wavelength observations for a substantial number 
of nearby SNe, irrespective of the type, galaxy host prop- 
erties, covering from as soon as possible after explosion 
until the nebular phase, is eagerly sought. 
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